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Abstract 

Feynman's laws of quantum dynamics are concisely stated, discussed in compari- 
son with other formulations of quantum mechanics and applied to selected problems 
in the physical òptics of photons and massive partides as well as flavour oscillations. 
The classical wave theory of light is derived from these laws for the case in which 
temporal variation of path amplitudes may be neglected, whereas specific experi- 
ments, sensitive to the temporal properties of path amplitudes, are suggested. The 
reflection coemcient of light from the surface of a transparent medium is found to 
be markedly different to that predicted by the classical Fresnel formula. Except for 
neutrino oscillations, good agreement is otherwise found with previous calculations 
of spatially dependent quantum interference effects. 
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1 Introduction 



In 1983, Richard Feynman was invited to give the first series of Alix J. Mautner 
Memorial Lectures at the University of Califòrnia, Los Angeles. These lectures were sub- 
sequently published as the book: 'QED The Strange Theory of Light and Matter' [1]. 
They are a masterpiece of popular science. The experiment Feynman chose to describe in 
detail was reflection from, and transmission of light through, glass sheets of various thick- 
nesses, at near-normal incidence . Closely analogous experiments were first performed 
and analysed by Isaac Newton [2] 1 . After a simplified analysis, the complete quantum 
mechanical calculation for the experiments is presented in Chapter 3 of [1]. This is done 
in terms of 'arrows' (vectors in the complex plane) that represent quantum amplitudes 
in one-to-one correspondence with paths where a photon is scattered from each atom of 
the sheet of glass. No equations are used. To the present writer's best knowledge, this is 
the only place, in Feynman's writings or elsewhere, where this calculation may be found. 
When Feynman had earlier discussed the reflection and transmission of light in 'The Feyn- 
man Lectures in Physics' [3] he had presented the Standard text-book analysis in terms of 
surface boundary conditions on the elèctric and magnètic fields of electromagnètic waves, 
leading to the Fresnel formulae for the reflection and transmission coeflicients. 

The principal aim of the present paper is to work out, in full mathematical detail, 
the QED predictions for some of the experiments discussed in [1]. It is found, perhaps 
surprisingly, that the reflection coeflicient of light at a vacuum/glass interface, calculated 
by the path amplitude method described in Chapter 3 of [1], is in marked disagreement 
with the Fresnel formula! 

The essential ideas whose applications were developed in [1] were already concisely 
stated [4] in 'The Feynman Lectures in Physics'. In describing the path amplitude inter- 
pretation of quantum mechanics in these lectures however [5] , Feynman prefered to discuss 
non-relativistic electrons rather than photons, in correspondence with his published work 
on the Path Integral approach [6, 7] that was restricted to non-relativistic quantum me- 
chanics. With the exception of the short passage cited above with the revelatory title 
'How it Works' light was always described in [3] in terms of electromagnètic waves, not 
photons. 

Feynman stated that the path integral formulation of quantum mechanics is not only 
the most fundamental, but also the easiest to understand. Why then has its impact, both 
on the teaching of physics 2 and on the presentation of quantum mechanical subjects in 
the research literature been so limited? This is true even for subjects where the path 
amplitude approach would seem to be the most natural one, such as physical òptics and 
quark flavour and neutrino oscillations in particle physics. The answer to this question 

1 These experiments actually studied interference efïects of light crossing a thin air film trapped between 
sphcrical and plane glass surfaces -the well-known 'Newton's Rings'. The quantum mechanical analysis 
of these experiments is essentially identical to that of those discussed by Feynman in [1]. The jacket of 
the first edition of this book showed a coloured photograph of thin film interference effects produced by 
a thin oil slick. 

2 An important exception to this is the work of E.F.Taylor and collaborators [8] who have stressed the 
insight into classical mechanics and Hamilton's Principle provided by Feynman's formulation of quantum 
mechanics. 
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is perhaps related to the fact that the situations in which the path amplitude formal- 
ism is both transparent and powerful -quantum predictions for correlated observations in 
space-time, which may be termed 'quantum dynamics', is far removed from the problem 
that was confronted and solved by the founders of quantum mechanics: Bohr, Heisen- 
berg, Schròdinger and Pauli. This problem was to describe the structure of àtoms and 
to calculate the rates of radiative transitions between different atòmic energy levels, a 
subject which could be termed 'quantum statics' since the detailed space-time structure 
of fundamental processes does not play any important role 3 . Indeed the latter subject is 
not discussed at all in the vast majority of text-books on quantum mechanics 4 . If future 
research physicists are never systematically taught Feynman's formulation of quantum 
mechanics it is not surprising that they remain forever ignorant of it. When some brief 
discussion of, say particle propagation in space-time, is included in Standard text-books, 
concepts such as spatial 'wave packets', relevant only to certain classes of physical prob- 
lems, are often introduced, in a prefectly general way, leading, in some cases, to erroneous 
conclusions. This point will be considered in more detail in the following section. 

Feyman's formulation of quantum mechanics is more fundamental and powerful than 
the earlier formulations of Heisenberg and Schròdinger because, as Feynman showed, 
starting from the path integral formula Schròdinger's equation and therefore all of 'quan- 
tum statics' is readily derived. However knowledge of Schròdinger's equation alone is, 
by far, not suficient to obtain correct predictions for space-time correlated events. This 
is demonstated by the different predictions for 'neutrino oscillations' to be found in the 
published literature and text books [11, 12, 13]. 

Before the work presented in the present paper, the present writer is aware of only 
two papers in which the Feynman path integral formalism is formally applied to optical 
problems. The first [14] concerns the refraction of light and the second [15] the òptics 
of beams of non-relativistic partides. This being said, it remains true that in the de- 
scription of almost any experiment in quantum òptics involving múltiple photon paths, 
path amplitudes are tacitly employed, usually called instead 'photon fields', 'wave func- 
tions' or 'state vectors'. The most detailed explicit application of the Feynman path 
integral method known to the present writer is in atòmic physics: the ' Photo detachement 
Microscope' [16]. 

In all of the optical experiments to be described in the following, as is the case in [1], 
the initial state is assumed to be a single excited atom which subsequently decays sponta- 
neously producing a single photon. This photon interacts with the experimental apparatus 
and later the same, or a different, single photon is detected. Also as in [1], for simplicity, 
only experiments where the effects of atòmic or photon polarisation are inessential, and 
may be neglected, are considered. 

3 Bohr made the following remark on this subject: 'For example the experiments regarding the excita- 
tion of spectra by electrònic impacts and by radiation are adequately accounted for on the assumption of 
discrete stationary states and individual transition processes. This is primarily due to the circumstance 
that in these qüestions no closcr description of the space-time behaviour of the processes is requircd.' 
[9]. 

4 This is not the case for [10], an introductory text book largely based on Feynman's space-time 
approach. As in Section 2 below, Feynman's rules for constructing probability amplitudes are explicitly 
statcd in this book. 
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The structure of the paper is as follows: In the following section Feynman's principies 
of quantum dynamics and their interpretation are reviewed. Thus the bàsic formulae used 
to perform the calculations presented in the subsequent sections are given. In Section 3 
the invariant propagator of a free on-shell particle is derived from the relativistic Feynman 
path integral, and the temporal propagator of a particle at rest from the time-dependent 
Schròdinger Equation. In Section 4 the classical wave theory of light in which only spatial 
'waves' are considered is derived from the Feynman path integral. The calculation, from 
first principies, of the refractive index of a uniform transparent medium is performed in 
Section 5. A new phenomenon 'refraction annulment' which may occur when the excited 
atom decays promptly after production, is predicted. The laws of reflection, refraction and 
linear propagation of light are derived in Section 6. A quantitative estimation is made of 
the spread of the photon paths around the classical trajectory of ray òptics corresponding 
to an extremum of the phase of the path amplitude. In Section 7 the calculation of 
the reflection coemcient of light at normal incidence at the interface of two transparent 
media, as qualitatively sketched in [1] is presented. In Section 8 time-dependent variation 
of the fringe visiblity in a Michelson Interferometer, using as a light source spontaneously 
decaying excited àtoms produced at a known time, is calculated. This type of prediction, 
testing essential temporal features of the path amplitude formalism, cannot be obtained 
using the classical wave theory of light. In Section 9 the path amplitude method is 
extended to the description of de Broglie matter waves and heavy quark and neutrino 
flavour oscillations. In each case a 'two probability amplitude' experiment is analysed in 
detail. These are: Young double slit experiments using photons or electrons, neutral kaon 
flavour oscillations in the processes n~p — > A(Kg, ]<£,), Ks,Kl — > 7r ± e T z/ and neutrino 
oscillations in the processes: n + — > un — > e~p. A summary and outlook are given 

in Section 10. 

As this paper is a long one, the reader is recommended to look first at the concluding 
section for an overview, returning later to any earlier sections that contain material of 
particular interest. 

2 Feynman's Formulation of the Laws of Quantum 
Dynamics and Their Interpretation 

In Chapter 3 of Volume 3 of the 'The Feynman Lectures in Physics' [3] Feynman 
set down four general principies for the quantum description of space-time processes. 
These are rules concerning the meaning and the method of construction of the Probability 
Amplitude that provides the quantum description of any space-time experiment. Any 
such experiment is defined by a fixed (prepared) initial state \i) and a fixed (measured) 
final state |/). The probability amplitude, Afi, is a complex number that is constructed 
by summing all Path Amplitudes that have the same initial state \i) and the same final 
state |/). The path amplitudes, PA fi, are, in turn, constructed by multiplying together a 
time-ordered sequence of Process Amplitudes. Feynman's four principies then define the 
physical meaning of the probability amplitude and specify how it is constructed: 

I The probability, P/j, to measure the state, |/), given the prepared state, is: 

Pfi = (2.1) 
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In general the probability amplitude depends on the space-time evolution of the state 
which is a function of the spatial geometry and other physical characteristics of the 
experimental apparatus as well as of those of various process amplitudes where partides 
may be scattered, destroyed or created. It therefore does not correspond simply to any 
observable property of either the initial or final state and so should not be confused with 
the wavefunction of any quantum system. 

The path amplitude is constructed according to the law of Sequential Factorisation. If 
the state \i) evolves in time via a series of unobserved intermediate states: \kj), j = l,n, 
the corresponding path amplitude: (f\k n , k 2 , k\\i) is given by a product of process 
amplitudes: (kj + i\Tj +1 \j) , where is a transition operator that evolves the state \kj) 
into \kj + i). These process amplitudes are either the space-time propagators of partides or 
the invariant amplitudes of particle decay or scattering 5 processes that are calculated in 
momentum space according to the Feynman rules of QED or the Standard Electroweak 
Model. Since the latter have no space-time dependence they appear only as constant 
multiplicative factors in the path amplitudes. It is understood that all unobserved particle 
states in the process amplitudes that are not members of the chain of sequential states 
\kj), j — 1, n that define the amplitude, are to be integrated over. Since the corresponding 
integrals are space-time independent that contribute only a multiplicative constant to -P/i 6 . 

II Sequential Factorisation: 

PA fi ={f\k n ,...,k 2 ,k 1 \i) = (f\T f \k n )...(k 2 \T 2 \k 1 )(k 1 \T 1 \i) (2.2) 

The third principle is Quantum Mechanical Superposition. If several different series of 
unobserved intermediate states are allowed by the experimental apparatus, the probability 
amplitude is given by the sum over all such series of intermediate states: 

III Quantum Mechanical Superposition: 

^ = E···EE</I*»»···»*2»*il < > ( 2 - 3 ) 

The fourth principle applies if the initial and final states are composite, constructed 
from the tensor products of sets of states: |^), \f^), \f^)--- but the path 

amplitudes linking the pairs \i^),\f^); \i^),\f^);... have no intermediate states in 
common: 

IV Composite Factorisation: 

^/W/(2)...i(')i(!)... =Y[PAfU) i u) (2.4) 

j 

5 The scattering process may be elàstic or inelastic. 

6 In the examples discussed in the present paper no account is taken of these integrations. The 
calculated space-time quantum interference effects are unaffected by this omission 



4 



An important special case of this principle occurs when a common prepared initial state 
gives rise to a number of separately measured final states; i.e. I^ 1 **), |^ 2 ^),... — > \i), in 
(2.4), but the paths linking \i) to I/*- 1 -*), \ f^),.. are distinct, with no common interme- 
diate states. This results, when quantum experiments containing identical partides are 
considered, in what is conventionally termed an 'entangled wavefunction'. Examples are 
the decays of ortho- and para-positronium: 

e + e~ (l 3 S 1 ) - 777 , e + e~ (l 1 ^) - 77 

In these and similar processes quantum statistics must be respected. In the case of 
the observation of both photons from the decay of para-positronium, two different path 
amplitudes are possible: 

PA(12) /(1)/(2)ï = (/( 1 )|T D | 7 i)(7i|T(P 1 )|7i)(7i|TpK) 

x(/ (2) |T z? |7 2 )(7 2 |T(P 2 )|7 2 )(7 2 |T P |ï) (2.5) 

PA(21) /(1)/(2)í = (/ (1) |T d |7 2 )(7 2 |T(P 1 )|7 2 )(7 2 |TpK) 

x(/ (2) |T d |7i)(7i|T(P 2 )|7i)(7i|TpK) (2.6) 

In these formulae, 71 and 71 are the two identical and indistinguishable decay photons, 
\i) the initial para-positronium state, and |/ 2 ) the final states of the photon detec- 
tion processes, occuring at the ends of the paths, P\ and P 2 , through the experimental 
apparatus. Tp and Tp are the transition operators of the positronium decay and photon 
detection processes respectively, while P(Pl) and T(P 2 ) are photon transition operators 
for the paths P± and P 2 respectively. Bose-Einstein statistics requires that the path am- 
plitude, PAj(i)j(2)j, for the experiment is symmetric under the exchange 71 <-> 72. Thus 

PA f{1)fWi = PA(12) fWfWi + PA(2ï) f{1)fWi (2.7) 

In an experiment where two identical fermions are both detected Fermi-Dirac statisics 
requires the path amplitude to be anti-symmetric under exchange of the two fermions. 
The similar path amplitudes PA{12)j{i)^2) i and PA{21)^i)^(2) i should then be subtracted 
rather than added as in (2.7), to yield the path amplitude PA^i)j( 2 ) i for the experiment. 

In practice, experimental set-ups do not usually correspond to the preparation and 
measurement of a unique intial state and a unique final state, but rather sets of such 
states: l = 1,2... (/); \f m ), m = 1,2... (P) are prepared or measured, respectively. 
The appropriate prediction to be compared with experiment is then: 



m l 



k n k 2 fel 



(2.8) 



The incoherent sums over / and m correspond to the density matrices of initial and final 
states respectively. (2.8) is the bàsic formula that is used for the calculations presented 
below. It is a simple iteration of a formula first given by Heisenberg in 1930 [17] which 
was later adopted by Feynman as the basis for his space-time formulation of quantum 
mechanics [6]. 

The principies I-IV above describe the rules for constructing probability amplitudes, 
but are devoid of any dynamical content. This is provided by a fifth principle, the Feynman 
Path Integral, which gives the physical prescription to calculate the probability amplitude: 
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V The Feynman Path Integral [18]: 



A fi = | (/|exp|i^^||i>nf[ü(4(í)) (2-9) 



' paths 



3 

I 

V j=l 



where 7 



and 



S= L(xi,X!,X2,X 2 , ...Xp,X p , ...)dt 

J ti 

PW = L lm(£ ^0)//...//^...%i^ 

where Xq, xi, ... denote successive positions along the path each separated by a small fixed 
time interval, e, and A is a normalisation constant. The latter will be determined for 
the case of a free, relativistic, on-shell particle in the following Section. The function L 
is the classical Lagrangian of the system whose quantum mechanical behaviour is to be 
described, and S the corresponding classical Action. 

In the following, (2.9) is not used per se to calculate the entire probability amplitude 
but rather the individual process amplitudes, that are combined according to Feynman's 
second principle, to obtain the complete path amplitude for the experiment under discus- 
sion. The most important of these process amplitudes is the Green function, or space-time 
propagator, for an on-shell photon or an excited atom. These are derived in the following 
section, the former from (2.9) and the latter, more conveniently, from the time-dependent 
Schròdinger equation, shown by Feynman [6, 7] to be equivalent, in the non-relativistic 
limit, to (2.9). In his book on quantum mechanics Dirac used, conversely, the Schròdinger 
equation to derive (2.9) [19]. 

The physical interpretation of Feynman's principies I-V will now be discussed. In 
typical treatments of 'The Interpretation of Quantum Mechanics' in the literature it is 
customary to attempt to reduce the problem to its bare essentials by avoiding explicit 
reference to any actual application. The physical meaning of Hilbert space vectors, or 
wavefunctions, is discussed in complete generality, the dynamical content being limited to 
the 'unitary evolution' of fixed energy solutions of the time-dependent Schròdinger equa- 
tion. The danger of such a procedure is that of over-simplification. Instead of interpreting 
how quantum mechanics describes actual experiments, what is done is to instead 'intèr- 
pret' a simple mathematical model, that does not necessarily reflect the actual complexity 
of a real experiment. To avoid this danger here, the discussion focuses specifically on the 
type of experiment treated in [1]: 

(i) An excited atom is created. 

(ii) The atom decays spontaneously, emitting a single photon. 

(iii) The photon interacts with the experimental apparatus. 

(iv) A photon (not necessarily the decay photon) is detected. 



7 Here x p = (x p , x p ,x p ) and the dot denotes a time derivative 
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First consider the simpler sequence (i), (ii) and (iv) only. The atom is created, it decays, 
and the photon propagates through free space to the detector. If the detector has a 
surface area, S, perpendicular to the space vector, f, drawn from the excited atom to the 
detector, and the atom is unpolarised, there is no need to use quantum mechanics (or 
classical mechanics either) to calculate the time-integrated probability that the photon 
will be detected. All that is needed is a knowledge three-dimensional spatial geometry. 
The answer is: eDS/(4nr 2 ), where ep, is the efficiency of the photon detector. In order 
to calculate the probabilty, SP(tD), that the photon will be detected in the time interval 
5tp> at Í£>, quantum mechanics is needed, but only to calculate the mean lifetime, Tg, of 
the excited source atom. Suppose that the excited atom is produced by the passage of a 
pulsed làser beam tuned to the frequency (Ei — Ef)/h where Ei and E f are the excited 
and ground state atòmic energy levels. The làser pulse is assumed to have a Gaussian 
form with variance Of. If the mean time of passage of the exciting pulse is to an d the 
production time of the excited atom is tp, the above-mentioned probability is: 



SP(t D ) 



7r xAirr 2 ) r<? J-00 
e D S 



TS C 



2a? 



dtpSt 



D 



47TTgr : 



■ exp 



1 



07 



--(t D -t ---^) 



St 



D 



(2.10) 



The only physical postulates necessary to derive this result are an exponental decay law 
for the excited atom and a constant velocity, c, for the decay photon. 



It will now be instructive to re-derive this result using Feynman's principies. The 
corresponding path amplitude is, using II: 



PA fi = A D <L>| 7 | 7 >< 7 |T|A><7|A|P)Ap 
The five process amplitudes have the following meaning: 



(2.11) 



Ap is the production amplitude for process (i) above. 

(-f\Ai\P} 8 is the process amplitude for the propagation of the excited atom from its 
production time tp to its decay time í 7 (see (3.20) below): 



(7|A|P)=exp 



(2.12) 



Ef and Ej are the 'pole center-of-mass energies' (see Section 3 below) of the atòmic 
states i and /, and rg is the mean lifetime of the state i. 

( 7 |T|Aj) is the invariant amplitude for the atòmic transition: Ai — > Af + 7. 

(D| 7 | 7 ) is the propagator of the photon from its production time í 7 to its detection 
time tp, (see (3.11) below): 

1 



1 r i 

(D| 7 | 7 ) = - exp --[(£ 7 (í D -í 7 ) -p 7 r] 



(2.13) 



The last member of (2.13) follows since the photon is a massless particle with con- 
stant velocity c: c = E^/p^ = r / {t D — í 7 ). 



8 The notation used for space-time propagators is that (f\P\i) is the amplitude for a particle (or other 
quantum object, such as an atom or molecule), P, initially at space-time point Xi, to be found at Xf. 
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A D is the detection amplitude for process (iv) above. 



Also required for the calculation is the density distribution, pp, for production of the 
excited atom and the density, po of detected final states: 



dp P 

PD 




(t P - *o) 

2a? 



dt f 



(2.14) 
(2.15) 



Using I, noting that í 7 = t D — r /c, and integrating over tp gives: 



5P(í D ) PA = | \PA fi \ 2 dp P St D 



2 Al 




TC T* J-oo 



exp 



1 7\ (ÍP-ÍO) 2 



"(ÍD-ÍP--) 
^5 C 



2a? 



dtpStD 



— exp 



1 . r a? . 

ÍD-ÍO ^ 

r 5 c 2t 5 



5í 



(2.16) 



where 



A^|A D | 2 |( 7 |T|^)| 2 |( 7 |T|^)| 2 |A P | 2 
Comparing with the (2.10) shows that, for consistency: 

e D S 



|^| 2 |< 7 |T|A>l1<7mA}|lAp| 2 = 



Atïts 



(2.17) 



This is perfectly reasonable since |(7|T|Aj 



T$ = H/ts, and \Ap>\ 2 — S so that the 



photon detection probability is proportional to the area of the detector. 

The above comparison has shown the consistency of the calculation based on Feyn- 
man's quantum mechanical principies with the 'common sense' result (2.10). This example 
will now be compared with the typical description of the free space propagation of parti- 
des to be found in text books on quantum mechanics, quantum field theory, or òptics. For 
example, Dirac writes down the wavefunction of a free particle in an energy-momentum 
eigenstate as [20]: 



ip(xyzt) = ao exp[- 



p ■ x 



a exp[--(Et -p x x - p y y - p z z)\ (2.18) 



where x = (ct,x) and p = (E/c,p) are the space-time and energy-momentum 4- vectors 
of the particle. Dirac then states that such a 'plane wave' wavefunction has no physical 
significance, since it predicts the particle to have the same probability to be at any point in 
space-time 9 . A similar description the photon in the above example would be a radial wave 
function of the form exp[— ip • x/h]/r which would predict that the photon has the same 
probability to be found at any distance from the source at any time. Such 'wavefunctions' 
are evidently meaningless and bear no relation to the actual physical situation correctly 
described by (2.10) or (2.16). Dirac indeed agrees that such a wavefunction is meaningless, 
but suggests that the solution to the problem is to replace the plane wave in (2.18) by 

9 Actually, since the wavefunction of (2.18) is not square-integrable, the normalisation constant ao 
vanishes, as does also the probability that the particle is within any finite volume of space-time 
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a spatial 'wave packet'. This is then claimed to give a correct quantum description of 
actual free partides. It is argued that since the plane wave in (2.18) is an eigenstate 
of E and p, that it is natural from the Heisenberg Uncertainty Relations: ApAx > h 
and AEAt > h that that both x and t should be completely uncertain. Since it is 
certainly true that actual partides are never produced in an eigenstate of E and p it 
is further argued that they must therefore be produced with spatial wavepackets that 
are coherent superpositions of different E, p eigenstates. Indeed the photons produced 
by the decay of the excited atom in the above example do have an energy distribution 
of non-vanishing width. Nevertheless, there is no associated spatial wavepacket. The 
absurdity of describing the photon in the above example by a wave packet of the type 
suggested by Dirac, becomes evident on examining the physical consequences of such a 
hypothesis. Suppose that t s has the typical value of 1CT 8 sec. The correponding spread, 
Ap, in the momentum of the photon is given correctly by the energy-time Uncertainty 
relation as cAp = AE ~ Ts = h/rs 10 . According the the momentum-space Uncertainty 
Relation invoked by Dirac, the width, Ax, of the associated spatial wavepacket should 
be: Ax ~ h/Ap ~ cr s ~ 3m. If such a wavepacket were actually produced in the 
decay of the excited atom in the above example, quantum mechanics would predict that 
it could be instant aneously detected at a distance of, say, lm from the production point, 
in contradiction to the predictions of (2.10) and (2.16) above. Using, say, a photo-diode 
to detect the photon, in the above example, its spatial position is easily measured with 
a precision of Ax ~ lmm. This implies that, at the moment of detection, the product 
AxAp for the photon is known with an accuracy 3 x 10 3 better than allowed by the 
Uncertainty Principle. The latter is evidently not applicable in this case. 

The point is that the relation AxAp > % applies either to simultaneous measurements 
of x and p, as in the examples orginally discussed by Heisenberg [21], or to the character- 
istic widths of the spatial and momentum distributions of a bound state wavefunction, or 
the wave packets which correctly describe a beam of massive partides such as electrons or 
neutrons. These distributions are related by a Fourier transform, from which the corre- 
sponding Uncertainty Relation is easily derived. In the example above, knowledge of the 
excited state of the atom provides prior knowledge of the momentum uncertainty of the 
photon as Ti/(cts) where Ts is the previously measured or theoretically calculated mean 
life of the excited atom. No measurement is required to know the photon momentum 
with this accuracy 11 

What is not taken into account in Dirac's discussion, and that of other text book 
authors is that, in quantum mechanics, real ('on-shell') partides propagate over macro- 
scopic distances in a classical manner. This fact is taken fully into account in the Feynman 
space-time formulation of quantum mechanics, and, as will be demonstrated in the exam- 
ples worked out later in the present paper, is crucial in deriving its predictions. As will 
be seen, many of the latter are identical to those of the classical wave theory of light, that 
are well verified experimentally. The conclusion is that the position of a particle that has 
already been produced must be defined in a classical manner in order to correctly calcu- 
late the correponding path amplitude. The energy-time Uncertainty Relation TsTs = h 

10 It is here assumed that the lower energy level of the atòmic transition that produces the photon is 
the ground state, with infinite lifctime 

11 Here incoherent line broadening effects due, for example, to the Doppler effect, are neglected, so the 
atom is supposed to be isolated and at low temperature. 
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of course plays an essential role in the derivation of the common sense equation (2.10): 
there is quantum uncertainty in the time of decay. Once a photon is produced, however, it 
propagates over macroscopic distances according to the classical law: Ar/Aí = c within 
each path amplitude. Adding the path amplitudes to construct the probability amplitude 
modifies, via quantum interference, this simple behaviour, and gives rise to the purely 
quantum effects of refraction, diffraction and interference that will be discussed below. 

It is possible to replace Dirac's wavefunction (2.18) by one which does give a meaning- 
ful space-time description of the photon produced in the decay of the excited atom. Since 
the detection process can be considered as a probe of the spatial position of the photon, 
the 'photon wavefunction' can be defined by omitting the amplitude A D in (2.11) and mul- 
tiplying by a suitable 5-function to describe the classical space-time propagation property 
of the photon. Thus the radial wavefunction of the photon is given by (2.11)-(2.13) as: 



The wavefunction of the photon depends not only on the spatial position and the time 
as specified by r and t, but also on the time at which it is produced í 7 , as well as the 
production time tp of the excited atom. Using the Born probabilty rule: 



and the density distributions (2.14) and (2.15), the wavefunction of (2.19) is easily seen 
to give, with a suitable choice of the normalistation constant a , the same result for SP(t) 



The reader may now be asking why, since the common sense formula (2.10) already 
describes the probability to observe the decay photon at any time, assuming only an 
exponential decay law and that the photon is a massless particle, what is the use of 
Feynman's formulation of quantum mechanics (that gives the same prediction via a more 
complicated calculation) or indeed of quantum mechanics in any formulation, except for 
calculating the mean lifetime of the excited atom? The answer lies in the description of 
space-time experiments where principle III (quantum mechanical superposition) plays a 
crucial role. This is the case in all the examples worked out in the present paper. The 
essential information provided by the path amplitude concerns not the spatial position of 
the particle (which is calculated classically) but the phase of each path amplitude at the 
instant the final state is measured. It is these phases which control the observed quantum 
interference phenomena. However, the particle does move in a classical manner along 
each path. As Feynman strongly emphasised, the difference between a classical world 
and the actual quantum one is that the different paths in one-to-one correspondence with 
different classical histories of the quantum system contribute together, in parallel in time, 
according to (2.3), to the probability amplitude that describes a single detection event. 
This effect is most simply and graphically illustrated by the Young double-slit interference 
experiment, to be analysed in detail in Section 9 below, that Feynman chose to exemplify 
this behaviour [22], but in reality not two, but an infinite number of paths each with 
its own distinct classical history, are combined, in parallel, when quantum mechanical 
superposition operates. 




<J(r - c(í - í 7 )) (2.19) 




as (2.10) or (2.16). 
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It is interesting to note that the actual physical description provided by Feynman's 
formulation of quantum mechanics contains important elements of several different so- 
called 'interpretations ' of quantum mechanics. These are: 



(a) 'Consistent Histories' 

(b) Everett's 'Many Worlds' interpretation 

(c) de Broglie's original 'Pilot Wave' theory 



Each of these interpretations is now briefly discussed in comparison with Feynman's for- 
mulation. 

'Consistent Histories' [23, 24, 25] are usually presented in the language of Hilbert space 
projection and trace operators. In the notation of the present paper, and considering the 
simplest non-trivial case, with a single intermediate quantum 'property' [24], consistent 
histories are given by summed path amplitudes similar to those in (2.3) above, except 
that the sum over intermediate states is not over the complete set allowed by a given 
experimental configurat ion, as is understood in (2.3), but is instead limited to an arbitary 
domain, D k . Probabilities are assigned these sets of paths according to a definition similar 
to Feynman's first principle, (2.1): 



where 



P fi (D k ) = \A fi (D k )\ 2 
pG(D) „ feM 

A fi (D k ) = ^(f\T f \k)(k\T k \i) 



(2.20) 
(2.21) 

(2.22) 



(2.20) corresponds to the history definition of Omnès [24] and Gell-Mann and Zurek [25], 
whereas (2.21) is the original definition due to Grifliths [23]. The amplitude Afi in the 
denominator of (2.21) is given by replacing the limited domain D k of (2.22) with the 
complete set of allowed intermediate states as understood in (2.3). Histories defined by 
different domains D' k) D'l are said to be 'consistent' if the probabilites calculated according 
to (2.20) or (2.21) add in a classical manner: 

P fl (D' k + Dl) = P fl (D' k ) + P fl (Dl) (2.23) 

Conditions for this to be true were first given in the projection operator and trace formal- 
ism by Griffiths [23]. They are equivalent to requiring that the interference terms between 
the corresponding amplitudes Af^D^) and Afi(D k ) vanish: 

Re [A fi {D' k )A fi {D'iy] = (2.24) 

This condition is satisfied if the difference between the phases of Afi(D' k ) and Afi(D k ) 
is an odd múltiple of 7r/2, i.e. if the corresponding vectors in the complex plane are 
orthogonal. The 'consistent' histories are then those for which the quantum interference 
effects resulting from superposition vanish. The 'histories' associated with the Feynman 
path amplitudes that describe actual experiments have, of course, no such restriction. 
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The relation of Feynman's formulation of quantum mechanics to Everett's 'Many 
Worlds' or 'Relative State' formulation, can be seen by comparing their descriptions of 
the Young double slit interference experiment. Suppose that, initially, photon detectors 
are placed immediately behind the two slits and an observer (with a lot of time at his 
disposal!) records the responses of the detectors when single excited àtoms are created 
nearby as a photon source. The observer will either note (mostly) that no photon is 
detected (0) or, more rarely, that a photon is detected either behind Slit 1 (1) or Slit 
2 (2). Then, as suggested by Everett, the results recorded by the observer in sucessive 
repetitions of the experiment can be stored as a string of numbers: e.g. 00010020001... . 
Another observer performing a similar series of experiments might see another sequence: 
00002000100... . For a finite number of repetitions of the experiment there are a count- 
able finite number of such sequences, each of which, Everett said correponds to a different 
possible 'history of the world', which will be hard-wired into the brain of the observer, 
or his recording device, at the end of his sequence of measurements. The actual history 
of any given observer is limited to just one such sequence, but all of the others describe 
'alternative worlds' each with the same weight as the one actually observed.. 

Consider now the simplest case (unlikely but not impossible) when the observer in 
the first world notes a photon behind Slit 1 in his first experiment whereas the observer 
in the second world notes a photon behind Slit 2. Suppose now that the detectors are 
removed from behind the slits and placed in a position where a photon passing through 
either slit may be detected by them. There is clearly now a correspondence between the 
path amplitude for the photon to pass through Slit 1 and the 'world' of the first observer, 
and the path amplitude for the photon to pass through Slit 2 and the 'world' of the 
second observer 12 . Thus in a certain sense Everett's two possible 'worlds' are simultane- 
ously present when the corresponding path amplitudes are added to give the probability 
amplitude for the Young double slit experiment. Since in any actual experiment an in- 
finite number of different path amplitudes contribute to the probability amplitude it is 
as though an infinite number of classically distinct 'worlds' contribute, in parallel, to the 
probability amplitude. This is just what Feynman called the only real deep mystery of 
quantum mechanics. 

In the early days of quantum mechanics, de Broglie 13 proposed an interpretation called, 
in French 'La Théorie de la Double Solution' but usually now refered to as the 'Pilot Wave' 
formulation of quantum mechanics [27]. The idea was that the space-time coordinates 
of partides should appear explicitly in the theory so as to describe as, in classical me- 
chanics, the 'real' positions of partides. There is also a purely mathematical, abstract, 
configuration-space wavefunction, ip, associated with the particle which is conjectured to 
'guide' the motion of the particle in such a way that that the statistical predictions of 
conventional quantum mechanics are recovered. This idea was taken up by Bohm and 
Hiley and developed into a fully-fledged deterministic alternative to non-relativistic quan- 
tum mechanics [28]. In this theory a 'quantum potential' is introduced that gives rise to 
classical forces that reproduce observed quantum behaviour. 

12 Photons passing through either slit do so whether or not there is a detector immediately behind it 
13 An extensive discussion of these ideas can be found in [26]. In a later development of the theory 
the particle coordinates were replaced by those of a local singularity in a 'physical' wave, u, that was 
conjectured to exist in addition to the usual abstract configuration-space wavefunction, tp. In this article 
de Broglie also imagined, like Feynman in [1], 'clocks' moving with the partides or waves and recording 
the phase information. 
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However, as seen above, in Feynman's formulation, the space-time coordinates of par- 
tides do, in any case, occur in an essential way in the equations of the theory. The physical 
interpretation of the photon wavefunction in (2.19) is of a photon moving like a classical 
particle according to the law Ar/Aí = c, but also carrying with it a certain phase, anal- 
ogous to de Broglie's 'Pilot Wave'. It is precisely this phase information, when combined 
with that of all other possible paths of the photon in any given experimental situation, 
that does in fact 'guide' the photon paths towards certain regions of space-time and away 
from others due to the effect of quantum interference. This process is exemplified in the 
calculations presented below in the present paper. 

Two closing remarks in this section devoted to general principies. The first is that in 
the classical limit, h — > 0, the sums over intermediate states in (2.3) reduce to a single 
term correponding to the classical path of the system. (2.1), (2. 2) and (2.3) then give: 

pClass = K/I^IO-^I^I^^I^K)! 2 

= K/|T>|^>| 2 ...KA5|T , 2|fc?>| 2 K*?|T 1 |i>| 2 

= PSn-P&PÏi ( 2 - 25 ) 

Here {k®)--- are states along the classical path that minimises the action S in (2.6). 
The last member of (2.25) is just a statement of Bayes Theorem 14 for combining condi- 
tional probabilities in classical statistics. 

Finally, a brief mention of the 'measurement problem' of quantum mechanics in re- 
lation to Feynman's formulation. It would seem, at first sight, that the traditional 
'problems' of the 'reduction of the wavepacket' superposition of Hilbert space vectors of 
macroscopic states [30], with the associated 'Schròdinger cat' paradox etc. are completely 
avoided in Feynman's formulation. The essential theoretical concept, the probability am- 
plitude is a complex number with a Lorentz invariant phase, which is completely and 
unambigously defined by Feynman's principies I-V and the experimental conditions that 
are to be described. On the other hand, Hilbert space state vectors are quite ambiguous, 
depending on the representation (Schròdinger , Heisenberg, Inter action,...) used to spec- 
ify them. Such a Hilbert space vector, unlike Feynman's probability amplitude, does not 
describe the results of actual experiments. ' Measurement s' in òptics do not produce pho- 
tons in quantum eigenstates, but rather destroy them. There is typically no 'wave packet' 
to contract, no superposition of Hilbert space vectors corresponding to different macro- 
scopic final states of the 'detector', one probability amplitude for a live cat, a completely 
different one for a dead one... . 

Feynman had very little respect for philosophers of science, and was too busy doing 
science to spend much time on such qüestions himself. This being so, the 'philosophy' of 
Feynman's formulation of quantum mechanics, unlike that of the Bohr, Heisenberg and 
Schròdinger one, still remains to be written. 



See, for example, [29]. 
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3 Relativistic Path Integral and the Propagator of a 
Free On- Shell Particle 



A suitable starting point for the discussion is an integral equation for the one-dimensional 
path integral [31]: 



f°° 1 

ip(x,t + e) = 7_ oo ^ ex P 



x-y 



if}{y,t)dy 



(3.1) 



where e denotes a small time interval, L is the relativistic Lagrangian of a free, on-shell, 
particle of pole mass, mp, and relativistic velocity, f3 = (x — y)/(ce), and the normalisation 
constant A 1 is to be determined by consistency in the limit e — > 0, y — > x. This formula 
is now generalised to three spatial dimensions by introducing the variable: ff = x — y, and 
changing the integration variable from y to ff: 



W,t + e) =///j- exp 



te / \r]\ 

h L \T 



ip(x — ff, t)d 3 r] 



(3.2) 



Choosing the origin of coordinates such that x = (0,0,0), the one-axis parallel to the 
velocity vector of the particle: ff = (rj, 0, 0) and making the replacement t + e — > t gives: 



^(0,0, 0,í) 



di] 2 / di] 3 



o A- 



exp 



xe 



ip(—r), 0, 0, t — e)c?77 



(3.3) 



where, without loss of generality, 77 has been chosen to be positive. The two and three 
axes are perpendicular to the direction of particle propagation. Since ip and L on the 
right side of (3.2) do not depend on i]2 and 773 the integral over these variables may be 
absorbed in the undetermined normalisation constant so that: 



}_ = I dï] 2 J dr] 3 
A~ A 3 

The relativistic free-particle classical Lagrangian is [32]: 

L[ - ) = -m P c\\\- 



(3.4) 



c 2 e 2 



—mpC 



(3.5) 



Making a Taylor expansion in 77 and e, of 0, 0, t — e), retaining only the zeroth order 

term and changing the integration variable from 77 to (5 gives: 



#),0,0,í) = J Q _ exp [--mpcy 1-/3 



V>(0, 0, 0, t)df3 



Hence 



ce 



dp + 0(ei],e 2 ) 



or, neglecting terms of 0(erj,e 2 ), 
Since 



A 

A = ce 



e = At = 7AT 



Ar 



(3.6) 

(3.7) 
(3.8) 
(3.9) 
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where t and r are the laboratory and proper times respectively of the propagating particle, 
then, using (3.8) and (3.9), as well as the relation ce = 77/ /3, (3.3) may be written as: 



f°° 3 

ip(0,0, 0,í) = / -exp 
Jo v 



tmpc 



Ar 



tp(—r], 0, 0, t — e)dr) 



= / ÜT(0, 0, 0, t; —77, 0, 0, í — Aí)^(-77,0,0,í- At)dr] (3.10) 
Jo 

where K(x 2 , t 2 ; x[, fa) is the covariant Feynman propagator of the particle between the 
space time points (x\,fa) and (x 2 ,fa). Exploiting the Lorentz invariance of the phase of 
the exponential in (3.10) the covariant propagator of a free, on-shell, particle may finally 
be written as: 



K(x 2 ,t 2 ·,x 1 ,t 1 ) = 



P 



\x 2 - Xi\ 

(3 



— exp 



trripc 



Ar 



\x 2 - Xi 



— exp 



- n (E(t 2 



fa) -p - (x 2 - Xi) 



(3.11) 



The propagator then has the form of a spherical wave, familiar from physical òptics. 
Although, in the derivation, the spatial interval 77 was assumed to be a 'small' quantity it 
should be noted that there is no length scale in the propagator of a free on-shell particle. 
In any case, Huygen's construction which, as will be demonstrated in the following section, 
is recovered in the path amplitude formalism, guarantees spherical wave propagation from 
the combination of, in principle, infinitely short spherical 'wavelets' [33]. Equation (3.11) 
is therefore of general validity. 

The mass parameter, mp, in (3.11) is, by definition, the 'pole mass' of the particle. In 
the case that the latter is unstable, with mean lifetime tm, (3.11) is modified to: 



K(x 2 ,t 2 ] X!,ti) 



P 



\x 2 - Xi\ 



■ exp 



ic 2 



(3.12) 



where the decay width T = Ti/tm has been introduced. The physical mass, m, of the 
particle is then smeared around the value mp according to a distribution of width ~ T/c 2 . 
This distribution (the momentum-energy space propagator of the unstable particle) is 
given by the Fourier transform with respect to Ar of the exponential factor in (3.12): 



K{E; E , T) = / exp 
Jo 



i , „ iF. A 

- 5 <* - y' At 



exp 



iEAr 



h ) 



) d(Ar) 



(3.13) 



where Eq = mpc 2 is the center-of-mass energy corresponding to the pole mass and E = 
mc 2 is the physical center-of-mass energy of the particle. The integral in (3.13) yields: 



K(E;E ,F) 



i(E - E ) - T/2 



(3.14) 



This propagator corresponds to the well-known Lorentzian 'natural' optical line shape 
in the case that the 'particle' is identified with the initial atòmic state in a spontaneous 
radiative decay process. For the decay of an unstable elementary particle it is the equally 
familiar Breit-Wigner amplitude. 
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It will be important, for the discussion of the damping of interference effects in Section 
8 below, to note that the imaginary part of the argument of the exponential in the space- 
time propagator (3.12) is a function of the pole mass, not the physical mass, of the 
particle. 

The dimension, [X -1 ], of the space-time propagator is a necessary consequence of its 
defmition in (3.10). Although the phase of this amplitude is Lorentz-invariant, the entire 
amplitude is not. The propagator may also be written as: 



fí(Ar,m P ,7) = exp 



impc 2 . 
' -Ar 



(3.15) 



where the appearence of the relativistic parameter 7 makes manifest the non-Lorentz- 
invariant character of the propagator. In the rest frame of the particle, (3.15) becomes: 



K(At, m P , 1) = exp 



im P c 2 AT 



(3.16) 



h 

which differs from Feynman's [34] asymptotic Lorentz-invariant space-time propagator: 



K asym (Ar,m P ) = mp A exp 
8-kicAt 



impc 2 . 
-Ar 



h 



(3.17) 



only by the factor mp and a normalisation constant. Neither (3.16) nor (3.17) is suit- 
able to describe the time evolution of the wavefunction of a particle at rest, due to the 
manifestly unphysical behaviour of the factor (cAr)" 1 originating in the spatial integral 
in (3.10) which is no longer defined for a particle at rest. In seems that the one dimen- 
sional propagator describing the time evolution of the state of a particle at rest can not 
be obtained as a simple limit of the four-dimensional space-time propagator. Instead, 
the propagator of a particle at rest may be simply and directly derived from the time- 
dependent Schròdinger Equation. For a particle at rest the non-relativistic limit necessary 
for the validity of the Schròdinger Equation is actually the same as the exact kinematical 
result, so that the propagator found in this way is both exact, and Lorentz-invariant. 
However the 'Schròdinger Equation' used here, though universally so-called in discussions 
of particle flavour oscillations, is in fact a relativistic version of this equation since the 
non-relativistic kinetic energy term of the conventional equation is replaced by the sum 
of this term and the rest energy of the particle. In the use of such an equation to describe 
atòmic transitions the rest mass term contributes an overall multiplicative phase factor in 
the solutions of the equation, without any physical consequences, since this phase factor 
cancels in all atòmic transition matrix elements. In physical òptics, as aleady mentioned 
above, the crucial dynamical element in the probability amplitude is the temporal propa- 
gator of the excited atom that constitutes the photon source, or, more precisely, the time 
dependence of the decay amplitude of such an atom. Although this formula may be found, 
for example, in Dirac's book on quantum mechanics [35], for clarity and completeness it 
is rederived here. 

Denoting the mass of the atòmic ground state by m and the pole excitation energy 
of the state by E® the Schròdinger Equation for this quasi-stationary state is: 



ih 



om 



^ = (m c 2 + EÏM) (3.18) 
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with the solution: 



hk(r)) = |^(0))exp 



t(m c 2 + E?) 

h 



(3.19) 



Hence the time evolution of the transition amplitude between the states \ipi) and \ipf) is 
give by the relation: 



^/(r)|T|^(r)) = (^ / (0)|T|^(0))exp 



% 



(3.20) 



where an imaginary part — iTi/2 has been added to Ef to take into account, as discussed 
above, the lifetime damping of the amplitude of the excited state. In path amplitude 
language, (3.20) may be interpreted as the product as the amplitude for the excited state 
to evolve from proper time zero to proper time r (the exponential factor) times the time 
independent decay amplitude of the state. The formula (3.20) is directly applicable to 
nuclear (3 transitions on making the relacement Ef — E® — > Ep where the latter quantity 
is the total energy release in the nuclear transition. As will be discussed in Section 9 
below, the temporal propagator of the neutrino source state is predicted by the Feynman 
path amplitude formalism to make an important contribution to the neutrino oscillation 
phase. The time dependence of the decay amplitude of a particle that is destroyed in the 
transition is given by the replacements Ef — > mpc 2 , Ej — > 0, in (3.20) where mp is the 
pole mass of the particle. That is, the 'excitation energy' of the particle is equivalent to 
its mass, while the mass-energy of the corresponding 'ground state' vanishes. 



4 Derivation of the Classical Wave Theory of Light 
from the Path Amplitude Formalism Of Quantum 
Mechanics 




Figure 1: Diffraction of a photon, of energy E { — Ef, produced in the spontaneous decay: 
Ai — > Af + 7 of an excited atom at S, hy a smaü hole, H, in an opaque screen. The photon 
is detected at D. The vectors f, r~[ lie in the plane of the figure 
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The path amplitude (f m \k n , ...k2, ki\it) defined in (2.2) above is specified by the initial 
and final states i\ and f m respectively, and by the set of unobserved intermediate states: 
k n , ...k2, k\. The latter may refer to different space-time positions of the same particle, in 
which case, the 'process amplitude' (k^T^kj) is the space-time propagator between kj : 
Xj,tj and ki : Xi,tj or, alternatively, to the transition amplitude of some quantum process 
in which a particle is scattered or new partides are created. Of particular importance for 
the discussion of physical òptics is the case of the creation of a photon in the spontaneous 
decay of an excited atom. In this case the first process amplitude (the rightmost one in 
(2.2)) describes the production of the excited atòmic state, Ai, from the ground state, 
Aq, at time t : 

(himu) = (Aí\Tp\Ao) (4.1) 

where Tp denotes an atòmic transition operator. The second process amplitude (/c 2 |T 2 |/ci) 
describes the space-time propagation of the excited atom from its production time ío, to 
the production time, í 7 of the photon. The third process amplitude: 

(k 3 \T 3 \k 2 ) = (iA f \T\A t ) (4.2) 

describes the decay of the excited atom into a photon, 7, of energy — Ei — Ef where 
Ei and Ef denote the energies of the initial and final atòmic states respectively. The 
appropriate formula to describe the product of process amplitudes (k 3 \T 3 \k 2 ) {k2\T 2 \ki) , is 
given by (3.20): 

(A; 3 |T 3 |A; 2 )(A; 2 |T 2 |A; 1 ) = {-yA f \T\Ai) exp 

The first factor on both sides of (4.3) is the (time-independent) decay amplitude of the 
excited atom in (4.2). The second is the temporal propagator of the excited atom, assumed 
to be at rest. 

The remaining product of process amplitudes {k n \T n \k n -i) ■■■{ki\T 4 \k 3 ) denote either 
free space propagation of photons according to (3.11) or, in the case of diffraction, refrac- 
tion or reflection, photon scattering amplitudes. The final (left-most) process amplitude 
describes the photon detection process: 

(fm\Tf\k n ) = {V m \T D \An ln ) (4.4) 

Here V m denotes the final state of the photon detection process 7„^4„ — > V m that is the 
same in all the path amplitudes with final state label m. This may be an activated atom 
in a photographic plate, a photoelectron, a conduction electron in a photo-diode etc. The 
detection process does not discriminate the state 7„ of the detected photon, and occurs 
at the same time, te, in all the path amplitudes with final state label m. 

To illustrate the application of the path amplitude formula (2.2), a simple case is 
considered where an exited atom is produced at S and the decay photon is detected at 
D after diffraction by a small hole, H, in an opaque screen (Fig. 1). The corresponding 
path amplitude is: 

AA = (P|T D | 7 >(Z)|7|tf>A i i//AS<tf|7|7> 

{-yAf\T\A i )('y\A\0){A i \Tp\Ao) (4.5) 



£°-tiy2)(t 7 -t ) 



(4.3) 
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The process amplitudes (7|^4j|0) , (ÍJI7I7) and (D\j\H) are, respectively, the space-time 
propagators of: the excited atom from its production time to to its decay time í 7 , the 
photon, before diffraction at the hole at time í# , and the photon after diffraction. The 
product of the three space-time independent process amplitudes , (Ai\Tp\Ao) , ( r yAf\T\Ai) 
and (T>\T D \~f) described above, is denoted below as A. The diffractive scattering ampli- 
tude, Adi/f, will be derived below. The parameter AS* is the area of the hole H, which 
is used as a weighting factor proportional to the number of space time paths crossing the 
hole. 

Since, for a massless photon, EAt — pAx = mAr = 0, there is no contribution to the 
phase of the path amplitude (4.5) from the photon propagators 15 . Using (4.2) and (4.3) 
the propagator of the excited atom is given as 16 



(P|A|0) = exp 

Noting now that 



l(£°-£°-;ry2)(í 7 -í ) 



(4.6) 



t D = í 7 + T -^± (4.7) 



the path amplitude of (4.5) may be written as : 



AA = —A díff ASe W \-Ul% — Ej — &(t D - 1 - ^) 
rr\ va l c 



(4.8) 



Here the factors 1/r, \jr\ originate (see (3.11)) from the photon propagators. 
Introducing the parameters: 



and 



Tic h 



where M P and Mp are the 'pole masses' of the initial and final atòmic states, the path 
amplitude may be written as the product of four factors. The first has no spatio-temporal 
dependence, the second depends only on r, the third only on r\ and the fourth only on 
the times to and t : 

AA = A- - ^dtfj^ _ e _ c ( m+p )( ÍD _ Í0 ) tt.U) 

T T\ 

In the description of diffraction or interference effects on the right of the screen in Fig. 
1, at a given detection time, to, and for sufliciently large vàlues of to — to, only the r\ 



15 This was noted by Feynman on P103 of [l].'Once a photon has been emitted therc is no further 
turning of the arrow as the photon goes from one point to another in space-time.' 

16 This is the physical origin of the 'stopwatch hancT discussed by Feynman in Fig 67 of [l]:'But when 
we construct a monochromatic source, we are making a device that has been carefully arranged so that 
the amplitude for a photon to be emitted at a certain time is easily caculated: it changes its angle at a 
constant speed like a stopwatch hand.' It is clear from (4.6) that the hand actually rotates in a clockwise 
direction, not an anticlockwise one as in Fig 67 of [1]. This change of sign has no effect on the predictions 
of the path amplitude calculations. 
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variation of AA is of importance. If the finite decay width of the excited atom may also 
be neglected, i.e., if p <C n, the T\ dependence of the path amplitude is that of a spherical 
wave field: 

U(n) = (4.12) 

n 

which satisfies the Helmholtz equation: 

V 2 U + k 2 U = (4.13) 

In this case, the description provided by the path amplitude formalism is mathematically 
identical to the classical wave theory of light. Diffraction is described by Kirchoff 's equa- 
tion [36, 37] and all the well-known interference phenomena of light, such as Huygen's con- 
struction [33], the Young double slit experiment and the Michelson Interferometer, may 
also be described purely spatially, like analogous wave phenomena of classical physics, 
such as transverse waves on a string or surface waves on a liquid subjected to the force of 
gravity. In coherent làser òptics where the effective lifetime of the photon source is very 
long 17 , /) « k, and so the classical wave theory of light is expected to be a very good 
approximation. It is therefore necessary to look elsewhere for the specifically temporal 
effects that distingish quantum mechanics from the classical wave theory. 

The diffractive scattering amplitude, Adiff, in (4.5) is derived by applying Green's The- 
orem, with suitable boundary conditions, to (4.13). This yields the well-known Fresnel- 
Kirchoff diffraction formula [36], and the expression: 

Adiff{oi, a±) = — — [cosa + cosai] (4-14) 

4:71 

where the angles a and a± between the normal to the screen and the vectors r and 
f[ are shown in Fig l. It is interesting to note that, in order to derive the Fresnel- 
Kirchoff diffraction formula in the classical wave theory of light, it is necessary to assume, 
in addition to Green's theorem, that the wave field U vanishes at large distances from 
the point under consideration. In the path amplitude description, this damping occurs 
automatically due to the factor: exp[— p(c(íd — to) — r — ri)] in (4.11). Since (see Fig.l) 
cito — to) > c(tn — í 7 ) = r + r±, the argument of the exponential factor is always zero or 
negative, producing strong damping of the amplitude for diffractive scattering at distances 
much larger than r or r\ since, in this case, c(t D — t ) > c(t D — í 7 ) r + r\. 

The limit of the formula (4.14) for forward diffraction, can also be be derived, directly, 
using the path amplitude method. This is done by considering rectilinear photon proga- 
gation (propagation of plane waves in the classical theory of light). Consider a photon 
produced by a distant excited atom at O and observed in a photon detector at D as shown 
in Fig. 2. The x-axis lies along the line joining the atom and the detector. Following 
Feynman [39] a plane surface perpendicular to the x-axis, distant x from the atom and x\ 
from the detector, where x >> xi, is imagined to be divided up into small areas, from each 
of which the photon trajectòries from the atom may be diffracted towards the detector. 
In particular, diffraction is considered from a ring of radius R\ and width ARi, centered 
on the x-axis. Using (4.8) above and noting that c(tn — to) — r — r\ = c(í 7 — ío) 

À 

AA = — exp [-ík (c(t D -t ) -r - r\) - pc(í 7 - í )] ^l (fó/ /2Tr J R 1 Ai? 1 (4.15) 



17 Defining a 'coherence length': L co h = ch/Ti — 1/(2/9), typical vàlues of L co h for a spontaenously 
decaying atom and a coherent làser source are 3m and 30km respectively [38] . 
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Figure 2: A photon produced by spontaneous decay of an atom at S is detected at D. The 
integral of the path amplitudes of photons diffracted or scattered in a ring of width Ai? x 
and radius R\ in a plane perpendicular to SD is constructed. In the case of diffraction, 
(4.16), the ring lies within an imaginary surface at distance x from the source. In the 
case of scattering (refraction) , (5.1), the ring lies in a thin sheet of transparent material 
of thickness 5. 



Since (see Fig. 2), r\ = x\ + R\, it follows that RiARx = riAri. Summing over all paths 
diffracted at the plane, the probability amplitude at the photon detector is: 

r rnax 

A = — exp [—ík (c(t D — t ) — r)] 2nAdiff / ex P [inri) — pc(í 7 — í )] dr\ (4.16) 

T Jx\ 

where an arbitary upper limit, r™ ax has been assumed for the T\ integral. Noting the 
inequality: 

c(í 7 -ío) > r + r 1 - x - xi (4.17) 

it can be seen that large vàlues of r\ will give an exponential suppression of the integrand 
in (4.16). In this case the r\ integral is readily evaluated by use of the Huygens-Fresnel 
Principle [33]. This states that the value of the integral is one half of the contribution 
due to the first half-period zone. With the change of variable; = kti, (4.16) may then 
be written as: 



A = 



À 

r 
À 



exp [~ík (c(t D - t ) - r) - pc(í 7 - í )] 



■kA, 



diff 



KXl+TT 



exp [iip] dip 



exp [—ík (c(t D -t ) -r - xi) - pc(í 7 - í )] 



KXl 

2-kíA, 



diff 



A 

X S D 



exp [—ík, (c(t D - t ) - x S d) - pc{tj - í )] 



K 

2mA, 



diff 



K 



(4.18) 



where, because of the small variation of t 1 over the first half-period zone, the factor 
exp [— pc(t 7 — to)] has been taken outside the integral and in the last line, the approxima- 
tion r + x\ ~ x + x\ = x$d where xsd is the source-detector distance, has been made in 
the exponential factor and r ~ xsd in the denominator. The path amplitude for photon 
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propagation from the excited atom to the detector is given by an equation analogous 
to (4.5), with the fifth and sixth process amplitudes (in temporal order) omitted. The 
corresponding probability amplitude is: 

À 

A = exp [-ÍK(c(t D - t ) - x S d) - pc(t 7 - t )} (4.19) 

XSD 

Consistency with (4.18) then requires: 

M, = ~ = («o) 

where À° is the de Broglie wavelength of the photon 18 , in agreement with (4.14) for 
a = «i = 0. . 



The Refractive Index of a Uniform Transparent 
Medium 



The calculation of the diffractive scattering amplitude presented above is easily adapted 
to calculate the refractive index of a uniform transparent medium in terms of the atòmic 
density, J\í, and the elàstic scattering amplitude, A sca t, of a photon from an atom of the 
material. Thus, in Fig. 2, the plane in which the photons are considered to be diffracted 
is replaced by a thin, uniform, sheet of transparent medium of thickness 5. Assuming 
an isotropic angular distribution for the scattered photons, the contribution to the path 
amplitude of photons scattering once from the àtoms in the sheet 19 is, by analogy with 
(4.16),: 

A(S) = — exp [-ÍK (c(t D -t ) -r - r x ) - pc(í 7 - í )] J\ÍA scat 5 \ \ dr\d$ x (5.1) 
r Jxx Jo 

where the integration variable has been changed from i? x to r 1 and where the azimuthal 
angle, 0i, around the x-axis, has been introduced. If the sheet is infmite in transverse 
extent the integral may be evaluated using the Huygens-Fresnel principle, as done for the 
case of diffraction above. As in all actual experiments the transparent material has finite 
transverse dimensions it will be found interesting, for the following discussion, to assume 
that the limit r™ ax is a function of a set of geometrical parameters, d^ that define the 
transverse extent of the sheet, as well as of 4>i- The examples of rectangular and circular 



18 Note that A° is defined, according to (4.9) and (4.20), in terms of the difference of the pole energies 
Ef and E® of the atòmic states Ai and Af, not their physical energies Ei and Ef. 

19 Note a similar calculation of the radiation field from a 'sheet of oscillating charges' in Section 30-7 
of [3] , in particular the contribution of the upper limit of the r integral. See also the analogous calculation 
of the refractive index of neutrons interacting in matter in Section 5.3.4 of [10]. 
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geometry are considered in Appendix A. In this case the equation analgous to (4.18) above 
is: 

A(S) = — exp [-ík (c(t D - t ) - r)) - pc(t 7 - t )] i ^'^ scat ^ 
r k 

x / {[expfi/cxi] - exp[ÍKr 7 [ iax (di, 0i)]}d0i 
Jo 

À r . , ,. . x \ \ ,. . „ 2iriJ\f Ascat 



exp [-i/í (c(Í£> - *o) - ^sd)) - pc(í 7 - í )] 5 (5.2) 

xsd k 

since the 0i average of exp[mr™ aa; ] vanishes due to ràpid phase variation. Adding the 
contribution from the path amplitudes of the unscattered photon from (4.19) gives: 



A' = A + A(6) = A 

where 



1 2mMA scat 8 



Aexp[ió(j)] (5.3) 



= 27rAA^, caf (5 = 2ttJV Ascat fx ^ ^ 

K K 

The last member of (5.3) is vàlid providing that 5 is chosen sufficiently small that 50 <C 1. 
In (5.4) / denotes the ratio S/(x 1 +x) ~ i.e. the fraction of the total distance between 
the exited source atom and the detector filled with the transparent material. 

It is interesting to note that the result (5.2) may also be obtained by restricting the 
T\ integral to the contribution of the first half-period zone. Thus the effect of random 
geometrical boundaries in a medium of finite spatial extent is the same as the use of the 
Huygens-Fresnel Principle in an infinite medium. This equivalence is used to simplify the 
calculation of the reflection coefficient at the surface of a transparent medium in Section 
7 below. 

Since the detection time tu and the production time t of the excited atom are fixed, 
as is the distance between the photon source and the detector, the change of phase of 
the path amplitude induced by paths that scatter once on the àtoms of the transparent 
material, implies that the apparent velocity, v(f), of the light, between the source and the 
detector, when the fraction / of the space between the source and the detector is filled 
with the transparent medium, is less than the speed of light in vacuum. 



In vacuum: 



On adding the transparent material: 



-Kc[t D -t --) (5.5) 



c 




x\ ( 1 5<fi. 

= -KC I t D - t - -{X + —J 



«(/). 

where, combining (5.4) and (5.6): 



(5.6) 



c 2TxMA scatt f ( v 

1 + 3 (5.7) 



v(f) K 
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Defining the refractive index of the transparent medium, n, as: 



n=- = -£-r (5.8) 
v v(l) 



and using the last menber of (4.20), gives the well-known formula [40]: 

1 



n = l + 7 ^(X° y ) 2 AÍA scatt (5.9) 



2tt 

It is interesting to note that in the calculation, based on the the geometrical configuration 
shown in Fig.2, leading to (5.6), the photon corresponding to each individual path am- 
plitude propagates at the vacuum speed c, but when the path amplitudes corresponding 
to single scattering processes are combined with those of unscattered photon the phase of 
the total amplitude is the same as that of a photon propagating in free space, but with a 
reduced velocity. This is clearly seen by inspection of (5.5) and (5.6) above. The apparent 
velocity, v(f), has a linear dependence on the filling fraction /: 

v(f) = {l-f)c + fv (5.10) 

This equation may be re-arranged to give: 

Av c — v(f) 



Av max c-v 



= f (5.11) 



which states that the change of effective velocity normalised to its maximum value is 
equal to /, i.e. it is proportional to the number of scattering processes contributing to 
the probability amplitude. 

It may be remarked that the above calculation of the refractive index has been carried 
out entirely in the language of partides moving in space-time. The de Broglie wavelength 
À° has been introduced in (5.9) only to show that the usual wave-mechanics result is 
recovered in the path amplitude approach. 

The calculation of the refractive index just presented was performed on the assumption 
that the the sheet of transparent medium is sufficiently thin that the approximation of 
the last member of (5.3) is vàlid, and that only a single scattering event needs to be 
considered. Combining (5.4) and (5.9) the phase shift produced by the sheet is: 

*=**£^ (5-12) 

7 

Thus, for the approximation 1 + ió(p ~ exp iS(f> to be vàlid, 5 is required to be a tiny 
fraction of A°, say less than 6 x 10~ 7 cm, for the case of a Sodium D-line where A° = 
5.9 x 10~ 5 cm. 



The calculation of the refractive index is now repeated, first considering the case of a 
slab of transparent medium of thickness much greater than A°, so that the approximation 
of the last member of (5.3) is no longer vàlid, but assuming that the value of to — is 
appreciably larger than the typical photon flight time between the source and the detector, 
and secondly, also considering a thick slab, but assuming that tr> — t Q may take any value, 
including ones near to the minimum value fixed by the source-detector distance and the 
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Figure 3: A photon produced by spontaneous decay of an atom at S, far to the left of 
the figure scatters at the points P 1} P 2 and P 3 inside a block of transparent material 
before being detected at D. The path SP1P2P3D lies in the plane of the figure. Various 
geometrical parameters used in (5.13) are defined. 

velocity of light in vacuum. In both cases p = is assumed, so that damping effects due 
to the finite lifetime of the source atom are neglected. 

In order to calculate the refraction of light in a block of transparent material of thick- 
ness, L, much greater than A®, it will be found necessary to sum over all possible con- 
figuations of múltiple scattering of the photon from the àtoms of the medium. In Fig.3 
the source atom is situated far to the left, so that the path of the photon from the atom 
to the front surface of the block may be taken to be a constant, equal to the distance, x, 
from the atom to the front face of the block. As in Fig.2, the photon detector, D, lies on 
the x-axis, which is perpendicular to the faces of the block. The distance from the center 
of the block to the detector is x. In Fig.3, the photon scatters at the points Pi, P2 and 
P 3 which lie in a plane containing the x-axis, QD, and are at distances from D, along this 
axis, of xi, x 2 and x 3 . The path lengths between the points P\ and P2, P2 and P 3 and P 3 
and D are r\, r 2 and r 3 . Other possible paths of the photon in the block are generated 
by rotating the point P\ about the axis Q 2 P 2 to give the series of points: P{,P", ,P{"... 
all distant r\ from P 2 and lying on a circle of radius R\ in the yz plane. By varying the 
value of R\ the integral is performed over the positions of all possible atòmic scattering 
events between Xi and x\ + Axi. Similarly, by considering the points P 2 ,P 2 ', ,P'-l' ■■■ an d 
P3, P3, , P'-l' ... lying on circles of radius R 2 and R 3 respectively, centered on the axes Q3P3 
and QD, and integrating over all vàlues of R 2 and R 3 , all possible atòmic scattering events 
in the intervals from rr 2 to x 2 + Ax 2 and £3 to x 3 + Ax% are summed. Including all paths 
with radii in the range P, to P, + APj, (i =1,2,3) the single scattering formula analgous 
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to (4.15), generalises, for the three-scattering case, to 20 : 



AAs = 7 — i - i l V exp 



- )H ( c(t D - t ) - (X + X + ^ - Xl) 



eXp[tKri] MA scatt RiAR 1 A ( f )1 Ax 1 



x 

ri 

x exp ^ r2 ] J^A scat tR2AR 2 A ( j )2 Ax 2 

x eXp[mra] MA acatt R 3 AR 3 A<j> 3 Ax 3 (5.13) 

Changing the integration variables to ri, r2 and r 3 , as in deriving (4.16) from (4.15) 
and performing the integrations, neglecting, as in (5.2), the contributions from the upper 
limits of the integrals, on performing the azimuthal integrations, due to their ràpid phase 
variation, gives: 

À 

AA 3 = = j exp 

[x + x + j - Xl) 



-ík (c{t D - t ) - (x + x + ^ 
x exp [ík (—xi + (xi - x 2 ) + (x 2 - x 3 ) + x 3 )] Ç^-^Ascatt ^ /^ Xl Ax 2 Ax 3 



~ exp [—ík (c(t D - t ) - x S d)] { ^^AscatA A Xl /± X2 /± X3 (5.14) 

XSD \ K / 

Where the approximation x + x + ^ — x 1 ~ x + x + ^ = x$d has been made in the 
denominator, where xsd is the source-detector distance. Integrating now over x±, x 2 and 
x 3 gives: 



A(T\ A r • (U +\ M (iïKNA sca ttX T +t A F +f A T +ï 

A(L) 3 = exp [-ik (c(t D - t ) - x SD )\{ / dxi / dx 2 / 



3 i·^a-L· 

dx 3 

xsd \ k / J X2 Jxz ^ x ~~ 

À r . . . . 1 {i2nAfA scaU L\ 3 

= exp[-iK{c{t D -t ) - x S d)\^7 [ (5.15) 

Xsd 3! V k / 

The limits of the x±, x 2 and x 3 integrals are determined by the inequalities (see Fig.2) 
xi > x 2 > x 3 . In (5.15) the general expression for the nested n-fold integral: 

rx+% rx+± rx+± jn 

l n — \ dxi / dx 2 ... / dx n = — (5.16) 

Jx 2 Jx 3 Jx-\ n\ 

is used. This formula is derived in Appendix B. For the case of n scatterings of the photon 
in the block, (5.15) generalises to: 

AfT\ À r • / /, . \ 1 (i2itNA scatt L\ n 

A(L) 3 = exp [-ik {c{t D - t ) - x SD )\ — (5.17) 

x S n n\ \ k J 



20 Note that in (5.13) only forward scattering processes are includcd. Back-scattering between two 
àtoms separated by distances > A° results results in ràpid variation of the phase of the path amplitude 
strongly supprcssing such contributions. The vanishing of the amplitude for backward scattering in a 
uniform transparent medium can also be understood as an effect of destructive interference between 
diffcrent path amplitudes, as described in Section 7 below. 
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Adding the contribution of the unscattered photon from (4.19) gives, for the fully inte- 
grated path amplitude (i.e. the probability amplitude) at the detector: 



A(L) 



tot 



X 



A + A(L) 1 + A(L) 2 + ... 
À 

exp [-ík (c(t D - t ) - x S d)\ 

XSD 

i27lAÍA sca ttL 1 

k 2! 



i27rAÍA scatt L^ 2 



K 



A 

x S p 

À 

X S D 



cxp 



ÍK [c(t D - t ) - X S D 



+ ... 



2nN A SC attL 



exp [-ík (c(t D - t ) - (x S d - L + nL))] 



(5.18) 



Where the refractive index, n, defined in (5.8), has been introduced. Defining, as in (5.4) 
above, / as the fraction of the total distance between the source atom and the photon 
detector occupied by the transparent medium: 

/ = — (5.19) 



X S D 



(5.18) may be written as: 

A(L)tot = exp -ík (c(t D - t ) - -^jj 



X S D 



where: 



v(f) 



(5.20) 



(5.21) 



l + (n-l)/ 

which may be compared with (5.6) and (5.7) above. The phase of the amplitude in (5.20) 
agrees with that in (5.6) above, derived by considering single scattering of the photon in 
a thin sheet. 



The calculation of the refractive index is now repeated taking properly into account the 
time interval At = tp — 1 between the time of production of the excited atom í and the 
time of photon detection tp. Referring to Fig.3 it can be seen that, for the configuration 
of excited atom and detector shown there, the minimum possible value of At is equal to 
the time of flight in vacuum, xsd/c, of the photon between the excited atom and the 
detector and correponds to the straight line path SQCD between the source, S, and the 
detector. It also corresponds to the case when the excited atom decays promptly after 
production. For later detection times, the photon may follow paths displaced from the 
axis SQCD but, considering, for example, the case of triple scattering, the upper limits 
on the r 1; r 2 and r 3 integrals in (5.14) will now be restricted by the value of As + x$d, 
which is the maximum path length allowed by the actual vàlues of to and í and the 
source-detector distance xsd'- 

As = c(t D - t ) - x SD = cAt - x S d (5.22) 

Note that a path of length As + xsd still corresponds to prompt decay. Decays occuring 
at later times, for the same value of to, have path lengths less than As + xsd- Inspection 
of Fig.3 shows that r±, r 2 and r 3 must satisfy the inequality: 

As > n + r 2 + r 3 - Xl (5.23) 
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Noting also that; 



T\ > Xi — X2 

r 2 > x 2 — x 3 
r 3 > x 3 

the following inequalities may be derived from (5.23): 

x\ — x 2 < r\ < As — r 2 — r 3 + X\ 
x 2 - x 3 < r 2 < As - r 3 + x 2 

x 3 < r 3 < As + £3 

In the present case, the definite integral derived from (5.13) above is: 

A(L) 3 = A 



(5.24) 
(5.25) 
(5.26) 



(5.27) 
(5.28) 
(5.29) 



f X ~r"2 f X+ ~2 f X ' 

exp [-ík (c(t D - t ) - x S D + ^1)] (27T M A scaU ) / dx 1 / dx 2 \ 

XSD JX2 JX3 Jx-i. 

As— T2— rz+xi rAs—r:j+X2 



dXr, 



X 



X 



I·ISS- 

Jxi —x 



exp[iK,ri\dri 



X\—X2 
AS+X3 

X3 



exp[ÍK,r 2 ]dr 2 



X2—XZ 



exp[inr 3 ]dr 3 



(5.30) 



Performing the integrals over r±, r 2 , r 3 , x±, x 2 and x 3 , as shown in Appendix C, yields 
the result: 



A{L\ 



À 

XSD 



exp [—ík (c(t D - t ) - x S d)\ 



X 



l_ ( i2itNA scatt L^ 3 
3! 



K 



{iA$y 



+ ïA$ - 1 exp[iA$] + 1 



where: 



A<£> = kAs 



(5.31) 



(5.32) 



The nested r 1; r 2 and r 3 integrals in (5.30) generalise in a straightforward manner to an 
arbitary number of photon scatterings. The expression for the probability amplitude in 
(5.20) above then replaced by: 



A(L) tot = 



A 



exp [—ík (c(t D - t ) - x S d)} 



xsd 

x |exp[ï/3L] + exp[iA$] 
(ipLf ( (ïA$) 2 



-iPL + (*A$ - 1) 



+ 
+ 



3! 



2! 



+ ïA$ - 1 



(ipL) 4 /(iA$) 3 (iA$f 



4! 



where 



3! 2! 

_ 2irNA acatt 



+ iA$ - 1 + 



(5.33) 



(5.34) 



For vanishingly small vàlues of A<£>, the term containing the large square brackets tends 
to the limiting value: 1 — exp [if3L] so that A(L) tot reduces to the vacuum path amplitude 
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between the source and the detector and the refractive effect of the block of transparent 
medium disappears. 



In order to discuss the probability amplitude for arbitary vàlues of A<3>, it is convenient 
to re-write the infinite series in the curly brackets of (5.33) separated out into real and 
imaginary parts. The algebraic manipulations necessary for this are discussed in Appendix 
D. The following expression is obtained: 



A(L) 



tot 



+ 



+ 



_A_ 

X S D 



exp \—ík (c(t D - t ) - x SD )] 



x h+pL[S + i(l-C)] + 



2! 



[[C-1 + A$S + i[S - A<£> + A$(l - C)}} 



3! 



A$ - S + A$(C - 1) + ^-^S 



+ i[C-l + 



2! 



4! 

+ i[A$ 
+ ...} 



1 - 



2! 
(A$) 2 



2! 



- C + A$(A$ - S) + 



2! 



(c-i) + 



(A$) F 
3! 



(A$) ; 



(A$) 2 

S + A$(C-1 + ^^) 



(A$) 2 (A$) 3 
1 ; -(A$ - 5) + 1 ; 



2! 



3! 



1-C 

(5 



where 



5 = sinA$ 



C = cos A$ 



This formula for the probability amplitude may be written in a more compact fashion by 
introducing a notation for truncated series expansions of trigonometric functions: 



S = 0, Si = A$, S 2 = A$ 

(A$) 2 



3! 



A $ -í^£ + ^£,...(5.36) 



3! 



5! 



Cq — 1, Ci 



2! 



x _ (A$£ + (A$) 4 



2! 



4! 



(5.37) 



The expression in the large curly brackets of (5.35) is denoted as F(A$, f3L) re f. It is 
the complex amplitude that multiplies the vacuum probability amplitude A to take into 
account the interaction of the photon with the àtoms of the refractive medium. Using the 
definitions in (5.36) and (5.37), the real and imaginary parts of F re f may be written as: 



ReK 



ref 



I - ILS + ( ^-(C + SS 1 -l)+ ( ^-(CS 1 - Sd) 



2! 



3! 



+ y^(i-cc 1 -ss 2 )+y^(sc 2 -cs 2 )+ 



4! 



5! 



+ (_ir ^r (1 " CCn - 1 ~ SSn) + { - ir w^y {SCn " CSn) 



+... 



(5.38) 
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ImF ref = p L {i-C) + ^-(S-CS 1 ) + ^-(CC 1 + SS 1 -ï) 



+ 



+ 



y^(cs 2 - Sd) + ^£(1 - CC 2 - SS 2 ) 



2n\ 



5! 

5C n _0 + 



\2n+l 



' (2n+l)P 



CC„ SSr, 



+ ... 



(5.39) 



(5.38) may also be written as: 



ReR 



(5L i 



ref 



(3L {(3Lf 



1 + /3LS+^(J3L)(C + SS 1 -1) + , >t 



+ 



+ 



+ 



f3L {/3Lf 



4 
+... 

2n 



3! 



CCi 



(CSt - SC,) 
CS 2 ) 



2n-l 



2n + 1 
+... 



2n- 1)! 
{2n)\ 



1 — CC n _i — S 1 ^) 



f-ir^^(.sa t -c , s' n )) 



(5.40) 



As will be discussed below, the vàlues of quantities (3L and A$ typically satisfy the con- 
dition j3L 3> A$. Inspection of successive terms in the series (5.39) and (5.40) shows that 
there is a one-to-one correspondence between terms with similar trigonometric coemcients 
(e.g. 1 - CCi - SSi in (5.39) and 1 - Cd - SS 2 in (5.40)). However, in (5.40), the 
(f3L) n /n\ terms of (5.39) are multiplied by the numerically large factor —f3L/n. The con- 
tribution of not-too-large vàlues of n to ReF re f is then expected to be much larger than 
the corresponding contribution to ImF re f. For large vàlues of n both series are expected 
to converge rapidly, since in this case, all the trigonometric factors tend to zero: 

CC n + SS n - 1 ~ C 2 + S 2 - 1 = 

SC n - CS n ~SC-CS = 
It is then to be expected that, in general: 

\ReF re f\ ^> \ImF re f\ 

which implies that the phase, of the amplitude F re f is close to zero, independently of 
the vàlues of f3L and A$. In this case the phase of the vacuum path amplitude remains 
unchanged by the interaction of the photon with the àtoms of the refractive medium and 
a phenomenon that may be called 'refraction annulment' is expected to occur. Refraction, 
i.e. propagation of light in the medium with an apparently reduced velocity c/n, only 
occurs when, as in (5.14), the contributions from the upper limits of the r integrals vanish 
on integration over the azimuthal angles due to their ràpid phase variation. The latter is 
induced by the geomètric boundaries of the refractive medium that determine the phase 
of the integrand. These are effectively randomly placed at the distance scales of the order 
of the wavelength of light. This point is further discussed in Appendix A. However, for 
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Figure 4: A simple experiment to search for refraction annulment. The side view is 
shown of a cylindrical glass block of refractive index n = 1.5, radius R and length L 
with an oblique face containing the points N and N'. Far to the left is an intense photon 
source consisting of excited sodium àtoms. Photons produced by spontaneous decay of 
the àtoms corresponding to the minimum detection time to for a given production time 
t of the excited atom, can reach the detector D only by paths near ON and for decays 
occuring promptly after production of the excited atom. In this case the phase change 
of the probability amplitude due to interactions with the àtoms of the block is strongly 
suppressed, the effective refractive index is very low, and so the photon can be transmitted 
through the oblique face and detected at D. For later detection times, photons may be 
produced at later decay times, and all photon paths between ON and PQN } are available. 
Normal refraction then occurs in the glass, and each photon is totally internally reüected 
at the oblique face (e.g. the path ONM). 
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small vàlues of As, the photon paths always remain away from the transverse boundaries 
of the refractive material so that the latter have no effect on the value of the integral over 
paths. This is the case in (5.33). 

To estimate the importance of the refraction annulment effect consider the simple 
experimental configuration sketched in Fig.4. A prism is constructed by cutting obliquely 
the face of a glass cylinder of length, L, and radius, R, at 45° to the axis of the cylinder. 
In Fig.4, the normal to the oblique plane face of the prism is in the plane of the diagram. 
The source atom is assumed to be at a large distance to the left along the x-axis, which 
coincides with that of the cylinder, so that all photons from the source have essentially the 
same time-of-flight to the front face of the cylinder. The photon detector, D, is situated 
on the x-axis at distance / from the front face of the prism. The largest possible value of 
As for all the r integrals to be limited by the production time of the excited atom, rather 
than by the geometrical boundaries of the prism, corresponds to a photon trajectory such 
as PQD. A simple geometrical calculation yields: 

A$ mra = — As max = — — + O(-) (5.41) 

If the exited atom is sodium in the upper state of the D-line transition, \® = \ D = 
5.9 x 1(T 5 cm. With R = 5cm and / = 200cm, As max = 625/im and A$ max = 6.66 x 10 3 
rad or 1060 periods. Assuming a refractive index for the sodium D-line of n — 1.5 and 
using (5.9) and (5.34) with L = 40cm gives 

(3L = ^ 1)L = ^ = 2.12 x 10 6 (5.42) 
Ad Ad 

Thus, in this case, (3L 3> A& max 3> 1, so that, according to the argument given above, the 
phase of F re f, and hence the corresponding refractive index, is expected to be very small. 
For a given geometrical configuration, the ratio of j3L to A<è max is a constant, independent 
of the dimensions of the apparatus. In order for the upper limits of the r integrals to 
be independent of the geometrical boundaries of the prism, the excited atom must decay 
earlier than a time As max /c = 2.2 x 10 _12 sec after production. Since the lifetime of the 
excited state producing the sodium D-lines is 5.4 x 10 _8 sec, a fraction ~ 4 x 10~ 4 of all 
D-line photons emitted into the sòlid angle of the detector is then expected to cross the 
prism following an unrefracted trajectory close to OND, as if they were propagating in 
free space. These are the àtoms that decay during the time interval: < t < As max /c 
after production. A simple experiment using photographic or photon counting techniques 
would be sufncient to confirm of invalidate this prediction. For larger vàlues of As (i.e. 
larger vàlues oi tu — to) such that the upper limits of the r integrals are determined by 
the geometrical boundaries of the prism, as in (5.14), the usual refraction phenomenon 
will occur in the medium and all the photons from the source will be internally reflected 
from the oblique face of the prism 21 so that the photons will follow trajectòries close to 
ONM in Fig.4. 



The critical angle of incidence for total internal reflection is 42° for n = 1.5. 
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6 The Laws of Reflection, Refraction and Rectilinear 
Propagat ion of Light, and Fermat's Principle 



refractive index n 2 




Figure 5: A photon produced by the source S is observed in a detector at D, which is 
movable in a plane distant d from the oblique face of a block of transparent material of 
refractive index n\ surrounded by a medium of refractive index . The position of D 
that corresponds to an extremum of the phase of the probability amplitude for photon 
detection at D is found to correspond to SneWs law of refraction. 



In Fig.5 is shown the trajectory, SPOD, of a photon from the decay of an excited atom, 
S, at a point sufficiently far to the left of the figure that all straight line paths from the 
atom may be considered parallel to the line PO. Over the path segment, PO, of length, /, 
the photon propagates in a transparent medium of refractive index n\. The surrounding 
space, back to the position of the excited atom, is assumed to be filled with a transparent 
medium of refractive index ni- The interface between the two media, passing through O, 
is a plane surface whose normal is at an angle 6j = n/2 — a to the segment PO (see Fig. 
6a), where, in the plane defined by PO and the normal, the angle between the projection 
of the interface and PO is a. The photon detector, D, is at a fixed distance, d, from the 
plane interface at a position defined by spherical polar coordinates r, 9, <fi. The cartesian 
coordinate axes are defined so that the z axis lies along the normal to the interface and 
the path segment PO is in the xz plane. 

By analogy with (5.18) the probability amplitude corresponding to detection of the 
photon at D at time t D is 22 : 



A A 

A = — — exp 

r\lr 



-hc t D li, ~~) + i$ ( r > 



(6.1) 



where 

$(r, l) = /t(n 2 r + nil) 



22 In this section it is assumed that c(tn — to) — xsd is sufíiciently large that there is no refraction 
annulmcnt effect. 
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The position of the detector D that maximises the photon detection probability is now 
determined by requiring that $(r, l) is stationary. As first pointed out by Dirac [20], 
and much emphasised later by Feynman [6], this condition ensures that neighbouring 
trajectòries have almost the same phase, giving an amplitude, by vector addition in the 
complex plane, that maximises the modulus of the probability amplitude Afi in (2.1) 
above. The phase is varied by parallel displacement of the trajectory to P'O', where the 
point O', in the xy plane is specified by the coordinates: x = i?cos0i, y = i?sin0i so 
that (see Figs.5 and 6) r — > r' and l — > V where 

r' 2 = (dtan6>cos0- i?cos0i) 2 + (rftan6»sin0 - iïsin^i) 2 + rf 2 (6.2) 
/' = / + R cos 0i cos a (6.3) 

The parameters 9, 0, R and <fii are now varied so that the r, / dependent part, of the 
phase in (6.1) is stationary. This requires that: 

dr' dl' dr' dl' n 

n2 dR +ni dR =n w, +n wr° (6 - 4) 

Differentiating (6.2) and (6.3): 

dr' d tan 9 , , , , R ,„ 

-7^ = ; — ( cos cos 0i + sm0sm0i) + — (6.5) 

oR r r 

dr' Rd tan 9 . N , >. 

(cos sm 0i — sm cos 0i) (6.6) 

di' 



dR 

di' 



cos 0i cosa (6.7) 
— i?sin0icosa (6.8) 



90i 

Substituting (6.5) and (6.7) into the first member of (6.4) gives: 

ri2(itan^(cos0cos0i + sin0sin0i) R 
h n 2 — + n\ cos 0i cosa = (6.9) 
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The solution to this equation with R = = is: 



( h ni cosa) cos0i = (6.10) 

Since rf/r = cos#, (6.10) gives: 

sin 6 = sin # = — cos et = — sin 9j (6.1 1) 

which is Snell's law of refraction. Substituting (6.6) and (6.8) into the second member of 
(6.4) and setting = gives: 

^.n 2 <itaní9 . . , . . 

R( ni cos a) sin 0i = (6-12) 

r 

which is verified both by the condition R = and (6.11). (6.10) and (6.12) show that the 
phase is stationary for an arbitary value of 0i provided that R = = 0. 

When a = n/2, so that the photon is incident normally on the plane interface, so 
that 9j = 0, (6.11) gives Qo = 0. The stationary phase condition then requires rectilinear 
propagation of the photon between S and D. 

The above calculation is easily adapted to the case of reflection at the interface, within 
the medium of refractive index n\, by the replacements in (6.12): 

— > li - R , n 2 —^n 1 

where Or denotes the angle of reflection relative to the inward normal at the interface. In 
this case (6.10) gives: 

rti(— síuOr + sinó 1 /) cos0i = (6.13) 
i.e., 0r — 0i, the law of reflection. 

The solution, (6.11), for the stationary phase with R = = and with determined 
by Snell's law, implies that the photon trajectory corresponding to the stationary phase 
of the path amplitude lies in the xz plane defined by the incident trajectory PO and the 
normal to the surface separating the refractive media. 

The sizes of the deviations to be expected from the classical trajectory corresponding 
to the stationary phase condition, as well as the spatial extent, within the media, of 
the regions where photon scattering processes contribute significantly to refraction or 
reflection will now be estimated using the path amplitude formalism. 

In the region of the stationary point, the variation of the phase is determined by the 
second order partial derivatives of the phase with respect to the trajectory parameters. 
For the angle of refraction o , Taylor's theorem gives: 



A$ = $ - $ síaí = i 



d^ídR\ 2 <9 2 $ íd^V 
\ ~d0j + 90? \ do) 



(A9 ) 2 (6.14) 

stat 
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From (6.1): 



<9 2 $ 
dR 2 


= K 


d 2 r' 
[ n2 dR 2 


d 2 l r 
+ Ul dR 2 _ 


(6.15) 


d 2 § 

«90? 


= K 


8 2 r' 

712 

d(f>( 


d 2 f 


(6.16) 



Taking the R derivative of (6.5): 



d 2 r' 



1 dr' 



R dr' 



dR 2 r' 2 dR 



— — d tan <9(cos cos fa + sin sin 00 + - - — — 



(6.17) 



Choosing now 0x = so that the photon trajectory is in the xz plane, and substituting 
the parameters at the stationary point: R = = in Eqns (6.5 ) and (6.17) gives: 



dR 
d 2 r' 



d tan 6q 



dR 2 

Combining (6.18) and (6.19): 

d 2 r' 



stat 



stat 



d dr^_ 
r dR 



stat 



tan#o H — 
r 



(6.18) 
(6.19) 



dR 2 



d 2 2n 1 

— — tan 2 6 + - 



stat 



r 



1-I-) tan 2 9 



o 



cos 2 9q 



The R derivative of (6.7) gives: 



d 2 r 

dR 2 







(6.15), (6.20) and (6.21) may be combined to obtain: 



d 2 



dR 2 



= ku 2 - 



cos 2 9o 



(6.20) 



(6.21) 



(6.22) 



stat 



Since cos 9 = d/r, 
d cos 9 



dR 



— sin 9q 



stat 



d9_ 
dR 



stat 



d dr 
'^dR 



stat 



d dr' 
'^dR 



sin 9q cos 9q 



stat 



That is, 



dR 
~dJ 



cos 9r 



(6.23) 



(6.24) 



stat — J 

Combining (6.14), (6.22) and (6.24) and noting that, since fa is chosen to be zero, the 
second term in square bracket in (6.14) vanishes, 



1<9 2 $ 
^dR 2 



stat 



'dR 
~d~9 



(A9 f 



Kn 2 r 



(A9 f 



(6.25) 



stat/ 
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With A$ = 7r, the spread of the angle do about its classical value is then estimated to 
be: 

A0 o = = \[^ (6.26) 

In a similar way, the displacements parallel to the x and y axes, Ax and Ay , corre- 
sponding to A$ = ir are found to be: 



Ayo = i/^-sec^o (6.27) 

V «2 



Ax = \ ^ (6.28) 

V n 2 

So, for the Sodium D-lines with À° = 5.9 x cm, refractive index of 1.5 and r = lm, 
and normal incidence, Axo = Ay = 6.3 x 10~ 4 cm. If the transverse dimensions of the 
interface are larger than this the photon paths are well represented by the classical light- 
ray corresponding to the stationary phase condition. For significantly smaller dimensions 
this is no longer the case and the photon paths are smeared over a much wider region. 
This is the diffraction-dominated domain of the classical wave theory of light. 

Introducing the effective velocity, v, of the photon in the refractive medium: v = c/n 
according to (5.6) the phase $(r, í) defined after (6.1) may be written as: 

$(r, i) = K(n 2 v 2 t e 2 ff + mvrfl 11 ) = Kc(t e 2 íf + tt ff ) = KcT efí (6.29) 

where 

f eff _ 1 f eff _ r 
Vl V 2 

and T e fï = t 2 + 1\ is the total effective propagation time of a photon from P to D in 
Fig.5. Since the propagation time from S to P' is a constant, independant of the position 
of P', a stationary value of $ correponds also to a stationary value of the total effective 
propagation time from the source to D. Thus the classical photon trajectory passing 
through D correponds to a stationary value of the photon's total effective propagation 
time, which is Fermat's Principle. 



7 The Reflection Coefficient of Light at Normal Inci- 
dence from a Plane Interface between Transparent 
Media 

A practical experimental set-up to measure the reflection coefficient of light is shown 
in Fig 7. A plane interface between vacuum and a medium of refractive index n lies in 
the yz plane of a Cartesian coordinate system. S is a light source (as previously a single 
excited atom produced at time í ) a large distance x = —x (x > 0) from the interface. 
In an actual experiment this could conveniently be replaced by a small source in the focal 



37 




X = -X r 




Figure 7: An apparatus to measure the reflection coefficient at normal incidence from a 
plane interface between vacuum and a block of uniform transparent material of refractive 
index n (see text). 



plane of a converging Iens or by a parallel làser beam. At (x,y,z) = (—Jo/2,0,0) is a 
half-silvered mirror (HSM) (or an equivalent optical component) MM, with the normal 
to its surface in the xy plane at an angle of 45° to the x-axis. A circular collimator CoCo 
restricts the photon paths to a cylindrical region around the x-axis. Two small photon 
detectors D\ and D 2 are situated at (— Jo/2, Jo/2, 0) and (— J /2, — 3Jo/2, 0) respectively. 
The virtual images of D\ and D 2 generated by reflection in the HSM, as viewed from O 
and S respectively, are denoted as D[ and D' 2 , and that generated by reflection of D\ in 
both the HSM and the surface of the transparent medium, as viewed from S, by D'{ . Since 
D' 2 and D" coincide (see Fig 7) it can be seen that, with this arrangement, equal sòlid 
angles are subtended at the detectors by the source S. Photons detected in D\ cross the 
HSM, are backscattered from the transparent medium and reflected at the back surface 
of the HSM (e.g. path SNPQRDi). Photons detected in D 2 are reflected from the HSM 
directly into D 2 (e.g. path SCD 2 ). 

It will be assumed that the decay width of the excited atom may be neglected (i.e. 
p = in (4.11)). The time dependence of the path amplitude for a given detection time 
t D is then ~ exp[— inc(t D — í )]. Assuming that the back surface of the HSM is half- 
silvered, as typically done in Michelson interferometers, which have a similar geometry, 
the probability amplitude for the photon to be observed in D 2 is 23 : 

A(D 2 ) = AoA D TH SM R H SM^[i{2(f) T HSM + 0§ 5M )] exp 

23 Note that a photon detected in D 2 crosses the HSM twice 



.3/tJo 



(7.1) 
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where AqAd denotes the probability amplitude for photon detection if D 2 were situated 
in the plane x = —/b/2, in the absence of the HSM, (including also the time depen- 
dent factor mentioned above), and Ad is the process amplitude for photon detection. 
Thsm ex P[^0/í5M] an d Rhsm ^WÍ^hsm} are ^e transmisson and reflection amplitudes of 
the HSM. 

The path amplitude, AA(Di) for the photon to be scattered at some point, P, in the 
interior of the transparent medium, and subsequently be detected in D\ is given by the 
product of seven process amplitudes: 

AA(D 1 ) = A D (D\>y\R)R H sMexp[í<i>% SM \{R\'y\P) 

yi{UA scat RARAxA4 l ){P\ 1 \C')T HSM eM^ T HSM\A ü (7.2) 

Here C is the intersection with the plane x = — Iq/2 of the photon path SN, and (R, (fi, x) 
are cylindrical coordinates specifying the position of P. Using (5.18) to write explicitly the 
triple product of photon propagators in (7.2), and integrating over all possible positions, 
P, of the scattering process gives 



ík[— + l — r + n(x + r)] 
2 



A{D 2 ) = I J I À! exp 
where 

Ài = AoA D R H sMexp[i(fi^ SM }T H sMeyip[i(j)HSM} 
In (7.3) the following spatial intervals have been introduced (see Fig 7): 

/ = PD[, r = PQ, x = NP = OL 

Noting the relation: 

lo l - r 



N AscatRdRdxdcfi 

/ [7 - ó) 



(7.4) 

which is a consequence of the similarity of the triangles D[QO and D[PL, the distance r 
may be eliminated from (7.3) to give: 



A{D 2 ) = 1 1 1 À! exp 



. MIq + ux) Iq 

lK[ n^x~ + 2 +nx] 



J\ÍA scat RdRdxd(fi 

/ (7 ' 5) 



As, according to the Huygens-Fresnel Principle, only vàlues of x of order À° contribute 
to the x integral in (7.3), it follows that, in the integrand of (7.3), x,r <^ Iq so that 
l 2 = R 2 + (Iq + x) 2 ~ R 2 + Iq is a good approximation. Thus RdR = Idl. Using this 
relation to eliminate R in favour of / in (7.5) and performing the / integration, in a similar 
way to the r x integration in (5.1), gives: 



J J Lhi\ In — IIX 



ík(Iq — nx)~ 



dxdcfi 



x < exp 



. Jmax((fi)(h + nx) lo 

lK \ ; l·" t: 

lo + x 2 



nx 



exp 



ík[— + 2nx] 



(7.6) 



The lower limit of the integral is l m i n = lo + x. On performing the (fi integration the 
contribution from the upper limit of the / integral vanishes due to the ràpid phase variation 
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resulting from irregularities of size À° or larger in the shape of the collimator Co that 
determines knaxi.4 1 )- 

Performing the x integral according to the Huygens-Fresnel Principle (the integral is 
equal to one half of the contribution of the first half-period in x), and neglecting x relative 
to / , except in the phases of the path amplitudes, gives the final result for A(D 2 ): 

Aín \ A (T^Ascat\ rMk, 

A(D 2 ) = - Aí {-^) exp[.— ] 



where, in the last line, (5.9) has been used. (7.7) gives the reflection coeficient at normal 
incidence, calculated using the Feynman path (FP) method, p^ p : 



P F R P 



n — 1 



2n 



(7.8) 



It is measured by observing the ratio of the counting rates of the similar detectors (assumed 
to be equally efflcient), D 1 and D 2 : 

RatejD,) _ \A(D X )\ 2 p F R p 

Rate(D 2 ) \A(D 2 )\* T' HSM { " > 

The modulus of the transition amplitude of the HSR, Thsm is readily measured by replac- 
ing the block of transparent material with a plane specular reflector of known reflectivity. 

In should be noticed that p R p differs markedly from the Fresnel prediction for the 
reflection coefficient at normal incidence: 



Fresnel 

Pr 



n — 1 



n + 1 



(7.10) 



For example, for a glass/vacuum interface, with n glass = 1.5, p^ p = 0.028, p^ esnel = 
0.040, a difference of 43%. A simple experiment similar to that sketched in Fig 7 could 
easily discriminate between these predictions. Another interesting difference with respect 
to the Fresnel formula is related to the minus sign on the RHS of (7.7). This sign, for the 
Fresnel prediction, is, in fact, ambiguous. Taking the limit as the angle of incidence tends 
to zero for photons linearly polarised perpendicular to the plane of incidence gives for the 
reflected amplitude (n — l)/(n+ï), whereas the similar limit for photons polarised in the 
plane of incidence is — (n — l)/(n + 1) [41]. Of course, at normal incidence, the terms 
'perpendicular' and 'parallel' polarisation become meaningless as the 'plane of incidence' 
is no longer defined. However, the actual phase shift for reflection at normal incidence 
is certainly measurable, and the path amplitude calculation, unlike the Fresnel formula, 
gives the definite prediction, n, for the phase shift. 

This phase shift could be measured by replacing the the detector D 2 by a suitably 
placed plane specular mirror and adding an optical attenuator so that the light reflected 
from this mirror and detected in D ± , and the light scattered by the transparent medium, 
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reflected by the HSR and detected in D\ have a similar intensity so as the maximise 
interference effects for a photon passing along either of the two possible paths. Such an 
experiment, using a suitably calibrated optical attenuator could also measure the reflection 
coefficient with a single photon detector. The so-modified apparatus is in fact a Michelson 
interferometer with one of the specular mirrors relaced by a vacuum/glass interface. The 
Feynman path amplitude description of a conventional Michelson interferometer is given 
in the next section. 

The result (7.7) for the probabilty amplitude is easily adapted to the case when the 
vacuum in Fig 7 is replaced by a uniform transparent medium of refractive index n\ while 
the refractive index of the original block of transparent medium is denoted by n 2 . Suitably 
modifying the photon propagators in (7.2), using (5.18) and calculating the contribution 
to the probability amplitude from photon paths scattered from the àtoms of the medium 
of refractive index ni, gives, instead of (7.7), the expression: 



which may be compared with the Fresnel formula, which gives, instead, for the quantity 
in the large curved brackets of (7.11): ±(n 2 — n\)/{ri2 + ni). The minus sign of the 
second term in the large curved brackets in the first line of (7.11) results from a reversal 
of the order of the limits in the x integration over the first half-period zone, since the 
àtoms of the medium of refractive index ri\ are situated at negative x. Setting n 2 = n , 
ni = 1 in (7.11) recovers the previous result (7.7). If rii > n 2 there is no phase shift on 
reflection in accordance with the well-known result of the classical wave theory of light. 
If rii = n 2 A(D 2 ) vanishes so there is no back-scattered light in a uniform refractive 
medium. According to (7.11) this may be interpreted as the result of perfect destructive 
interference between the path amplitudes of photons back-scattered scattered from àtoms 
with positive and negative x. 

Interference phenomena closely related to those occuring in 'Newton's Rings' [2] that 
are extensively discussed in [1] are simply analysed by suitably modifying the x integration 
limits in (7.6) above. Replacing the block of transparent medium by a thin sheet of 
thickness À°(l + 2p)/An, where p is a positive integer, gives a reflection coefficient four 
times larger than in (7.8). If the sheet has a thickness À°(l +p)/2n, p^ p vanishes due to 
perfect destructive interference, in this case, of the path amplitudes corresponding to the 
scattering of the photon from each atom of the medium comprising the sheet. 



8 Spatio- Temporal Interference Effects in the Michel- 
son Interferometer 



A schematic layout of a Michelson interferometer [42] is shown in Fig 8. The photon 
source, S, a single excited atom, is produced at a known time, to- For concretness, as 
discussed in Section 2 above, this can be taken to be an atom of, say, sodium vapour, 




(7.11) 
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M, 



L + d 




i/* 

D 



Figure 8: The Michelson interferometer. A photon from the decay of an excited atom at 
S, produced at time Íq, foüows the paths SOMiOD or SOM 2 OD in alternative histories 
of the decay process and subsequent space-time propagation of the photon, before being 
detected at D at time t D . PP is a glass plate half-silvered on its back surface. S[ is 
the virtual image of S after reüection fírst in PP then in the plane mirror M-y. S' 2 is 
the virtual image of S after reüection first in the plane mirror M 2 then in PP. The path 
diffrence 2d is chosen to be > crg where t$ is the mean lifetime of the excited source atom. 
Note that the atom must decay at different times in its alternative histories in order that 
the photon arrives at the detector at the fixed time to- It is just this time difference, 
determining the phase advance of the propagator of the excited atom, that determines 
the size of the quantum interference term generated by the amplitudes corresponding to 
the two alternative paths. 
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resonantly excited by a tuned, pulsed, làser beam. The pulse width should be at most 
a few nanoseconds and the time of passage should be known with a precision of lns or 
better. 

The photon produced by decay of the excited atom may follow the paths SOMiOD 
(reflection in by the mirror M ± ) or SOM 2 OD (reflection in by the mirror M 2 ) before being 
detected, at time to, in the photon detector at D. The glass plate PP situated at the origin 
of a Cartesian coordinate system, contains the z-axis of the latter and is inclined at 45° to 
the x- and y-axes. The back face of the plate is partially silvered to increase the moduli of 
the path amplitudes associated with the paths SOMiOD and SOM 2 OD. In a practical 
interferometer, focussing lenses would be installed in the paths SO and OD to render 
the photon trajectòries in the interferometer parallel to the x- and y-axes and increase 
the efficiency of photon detection. They would contribute a constant multiplicative factor 
to the interfering path amplitudes and so have no effect on the interference phenomena 
discussed here. 

The simplest way to understand and analyse the operation of the interferometer is to 
consider the virtual images S[ and S' 2 of the source, as reflected in the partially silvered 
plate and the mirrors M\ and M 2 , respectively, and viewed from the position of the photon 
detector [42] (see Fig 8). The following length intervals are defined: 

SO = OM 2 = OD = L, OM 1 = L + d 

It then follows that: 

S[D = Li = AL + 2d, S' 2 D = L 2 = 4L 

The probability amplitudes, .4.(1), -4(2) for photons reflected at Mi,M 2 and detected by 
D at time tn are given by (5.18) as: 

-4(1) = ^ex P H(/ ÍC --^)(í D -í --)] (8.1) 
Li 2r s c 

A(2) = ^expH(«c--^-)(í D -ío-— )] (8.2) 
L 2 2t s c 

where the the space-time independent amplitudes A± and A 2 include the production and 
detection process amplitudes as well as the amplitudes describing all transmission and 
reflection processes in the arms of the interferometer. In (8.1) and (8.2) the mean life 
T s — ^/r,s of the source atom has been introduced as in (3.12). For the Sodium D-lines 
Ts — lOns. This value will be taken for the quantitative predictions presented below. 
Using (2.1) and (2.3), the probability, P(t D < t 1 ^^) to detect the photon in D during the 
time interval t < t D < tp ax is given by the following expressions: 

(i) For: t™ x <t + L 2 /c 

P(t D < t^ ax ) = (8.3) 

(ii) For: t + L 2 /c < t% ax < t + Ljc 
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p(t D <t^)= r \a{2) \ 2 dt D 
jío+í/2/c 

(iii) For: í + L x j~ 



r S \À 2 \ 2 

a 



c < t^ ax 



l_exp[-l(í---í -— ) 

T S C 



(8.4) 



^max 

P(t D <t™ x ) = D \A(1) + A(2)\ 2 dt D 

J t rnzn 

£max -f-max 

= / ° \A(2)\ 2 dt D + / ° (\A(1)\ 2 + 2Re[A(l)A(2)*])dt D 

Jt +L 2 /c iío+Li/c 



= r s 



+ 



X 



X 



1 



\À 2 \ 2 

exp[- 



exp[ {t^ x 



[ -)\ 



max 
D 



c 



-4i| 2 , ol^ill^l 

+ 2 — ^—^ — exp 



L 1 L\L 2 
cos{k(L 1 - L 2 ) H 



2cr s 



'12 



}]} 



15) 



where 



012 = phase(^li) — phase(v4.2) 

For condition (i) the photon can arrive at the detector by neither path, so its detection 
probability vanishes. For condition (ii) the photon can arrive at the detector only via 
the p&th.SOM 2 OD so that there is no interference phenomenon. For condition (iii) the 
photon may arrive at the detector via either path so interference is possible. With the 
aid of optical compensators in the arms OM 1 , OM 2l A\ and A\ may be chosen so that: 



A ± = A 2_ =K 
L\ L 2 

In this case (also setting t = 0) (8.5) simplifies to : 

P(t D <t™ ax ) = t s K 2 {2- f(t™ x )(l + e^) 

+ 2(e~^ - f(t^ ax )e^) cos(2/trf + 12 )} 

where 

f(t) = exp 



(8.6) 



(8.7) 



T S C 



Neglecting small corrections due to the exponential terms in (8.7) the maximum and 
minimum vàlues of P{to < t"^i ax ) are: 



2d 

Pma X (min) T S K 2 {2 - f(t™ X ) (1 + e~ ) 

+ (-) 2(e"^ -f(t™ ax )e^)} 
and the 'fringe visibility' V(t D < t r £ ax ) is: 

Pmvc ~ Pmin 2(e~^ - /(í^)e^) 



V(t D < Q ax ) = 



+ Pn 



2- f(t 1 g ax ){l + e"s) 



1.8) 



(8.9) 
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Curves of V(t D < t r ^ ax ) as a function of tp ax are presented in Fig 9 for: ts = lOns, 
L = 50cm and d = 12.5, 25 and 50cm. Such interference effects as a function of t^ ax — 1 , 
are readily observable. It is suficient to adjust the interferometer for maximum destructive 
or constructive interference near some fixed value d = do and measure V, by small 
variations of d around d , for different vàlues of tp ax . If N max (N min ) are the numbers 
of photons recorded with t D < when the interferometer is adjusted for maximum 
constructive (destructive) interference then: 

MÍ+ ^ +max\ N max — N m i n 
V{t D <t D )exp = J7 

which may be compared with the prediction of (8.9). Fitting this prediction to the data 
determines the the only parameter, t$, that is not fixed by the geometry of the experiment. 



1 




Figure 9: Fringe visibilty V in the Michelson interferometer of Fig.8 as a function oft^ ax 
the maximum observation time of the photon at D. The excited source atom is produced 
at time t = 0. L = 50cm and t s = lOns. The curves A,B and C correspond to d = 12.5, 25 
and 50cm. 



Choosing a value of t™ x cr s results in f(tp ax ) ~ 0. The time-integrated fringe 
visibility as a function of d, V^i^d), is then qiven by (8.9) as 

^oo(rf) = V(t D < oo) = exp (8- 11 ) 

Thus the time integrated visiblity is predicted to decrease exponentially with the path 
difference 2d. Exactly the same behaviour is predicted by the classical wave theory 
of light [43]. In this calculation, similar to that of Michelson [44], The visibility was 
calculated using the equation obtained by making the substitutions k — > 2ticu 1 and 
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j-^rnax^ ^ CTs - m (g # 7) a bove, and by weighting the cosine interference term with a 
Lorentzian distribution in the freqency, v 1 , of the photon. In the Feynman path am- 
plitude calculation k is defined by the pole masses of the initial and final states (see 
(4.9) above), which, unlike the photon frequency, do not vary on an event-by-event basis. 
Thus, from the point-of-view of the Feynman path amplitude calculation, the classical 
wave calculation is incorrect, even though it predicts the same result. 

In the above treatment, two potentially important physical effects have been neglected. 

The first effect is 'pressure broadening' of the atòmic linewidth. This results in an 
observed value of Ts less then that, Tg at , corresponding to the 'natural' linewidth of an 
isolated, freely decaying, atom: 



the source of excited àtoms tends to zero. The physical origin of this effect is easily 
understood in the path amplitude language. In order for the source atom to retain the 
coherent phase relationship between different path amplitudes, following from (4.3), that 
is assumed in (8.1) and (8.2) above, the atom must remain in the excited state sufficiently 
long that unhindered spontaneous decay can occur. If this is not the case, for example, if 
the atom is dexcited, or excited into a different state by inter-atomic collisions, after the 
time corresponding to the 'earlier' path amplitude but before the time corresponding to 
the 'later' one, no interference will be possible. In the case that all inter-atomic collisions 
destroy the excited state then the parameter rp is simply related to the mean time between 
such collisions. There is then competition between decay and inter-atomic collisions for 
destruction of the excited state. On the assumption that these are independent processes: 



Which is equivalent to (8.12) above. In fact, the actual situation is much more compli- 
cated, since the inter-atomic collisions do not always result in destruction of the excited 
atom and may modify thee energy level or the decay probability of the excited atom. 
Many attempts were made during the '30s and '40s of the last century to calculate Tp 
from phenomenological atòmic models [45, 46]. 

The second effect is the possible correction due to motion of the source atom, which 
has previously been assumed, throughout the present paper, to be at rest. An analogous 
calculation has recently been performed [12] by the present author for the related 'neutrino 
oscillation' problem. It is adapted to the present case in Appendix E below. The result is 
markedly different from 'Doppler effect' formulae found in the previous literature and text 
books. There is no damping effect at any order in v/c, only a phase shift. The classical 
wave theory of light predicts a Gaussian dependence of the fringe visiblity on the path 
difference due to the first order Doppler Effect (DE): 



This formula, due to Lord Rayleigh, is obtained by weighting the cosine interference term 
in the visibility function by a Gaussian distribution of z/ 7 derived, using the Doppler 




(8.12) 
pressure of 




T s = T n s at + T P 



(8.13) 




(8.14) 
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Transition 


A? (A) 


A exp (cm) 


A DE (cm) 




t s (ns) 


rg at (ns) 


t p (ns) 


H r 3p-2s 


6563 


19.0 


14.3 


225 


0.46 


5.4 


0.50 


H 6 4p-2s 


4861 


8.5 


10.6 


516 


0.204 


12.4 


0.207 


Na D 3p-3s 


5893 


80.0 


65.0 


225 


1.92 


12.4 


2.98 



Table 1: Measurements from [44] of path lengths A exp yielding a fringe visiblity of 50% in 
comparison with theoretical expectations A nat derived using (8.11) and A DE from (8.14). 
The corresponding vàlues of the lifetime prameters ts and Tg at are also shown as well as 
the pressure-broadening parameter Tp derived from (8.12). 



effect, from a Maxwellian distribution of source-atom velocities. As in the classical wave 
derivation of (8.11) this requires the substitution /t — > Itïcv^ in (8.7) above, which is 
incorrect for the Feynman path amplitude calculation. 

The prediction of (8.14) was compared by Michelson in [44] to a number of ex- 
perimental measurements of fringe visiblity. The results obtained for some transitions in 
Hydrogen and Sodium are summarised in Table 1. The experimental observable, A, is the 
path difference 2d in the Michelson interferometer for which V 00 (^) = 0.5. Also shown in 
Table 1 are the measured vàlues of t s = A/(2c/n2) derived from (8.11), A DE from (8.14), 
theoretical vàlues of Tg at 24 and vàlues of the pressure-broadening lifetime Tp calculated 
from Ts and Tg at using (8.12). At the time that Michelson performed the measurements 
shown in Table 1, Quantum Mechanics had yet to be invented so no theoretical vàlues 
of Tg at were available. The vàlues of t p shown in Table 1 show that pressure-broadening 
effects are very important for the Hydrogen lines, but less so for the Sodium D-lines. 
Michelson had observed that line widths increase with increasing pressure, but compared 
his measurements only with the Doppler formula (8.14). For this comparison, in the case 
of the Hydrogen lines, Michelson chose a temperature of 485°K in (8.14) corrsponding to a 

R.M.S. velocity, \J kT/M, of 2.0 x 10 3 m/sec. Equation (8.14) was used to predict vàlues 
of A/À° for different atòmic transitions, which were compared with the corresponding 
experimental quantities. Rough agreement was claimed although some positive devia- 
tions of up to 30% and negative deviations up to 40% were observed in some cases. No 
uncertainties were quoted on the experimental measurements. Given the importance of 
pressure broadening effects for the measurements presented in [44], and the difficulty to 
distinguish between the former effect and that due to source motion, these measurements 
provide no evidence for or against the possiblity of a vanishingly small source-motion 
correction, as predicted by the path amplitude calculation. 

Precise measurement of the temperature dependence of the line widths of excited àtoms 
in a gaseous source then constitutes another stringent experimental test of path ampli- 
tude predictions. In practice it was very difficult in the past to disentangle conjectured 
'Doppler effect' and actual 'pressure broadening' effects in the data. Both are expected to 
result in larger linewidths at higher temperatures for a constant volume gaseous source. 
One possible approach is to measure Ts in a Michelson interferometer, using a constant 
temperature gaseous source, as a function of pressure at low pressure. Extrapolating to 

24 The vàlues of Tg at for the hydrogen lines are taken from P136 of [47]. The value for the sodium 
D-lines is calculated using the formula: l/rg at = Airr e v mn 2f / '(3c) [48] where r e is the classical electron 
radius and / = 0.9755. 
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zero pressure trien gives a line width with contributions only from the natural width and 
possible source motion effects, as rp is infinite at zero pressure. The prediction of the 
Feynman path amplitude calculation is that the same extrapolated value of t$ should be 
obtained, by this procedure, independently of the temperature of the source. 



9 Spatially Dependent Interference Effects in Quan- 
tum Systems with Two Probability Amplitudes. 

In this section the Feynman path amplitude description will be applied to some analo- 
gous physical systems: the Young double-slit experiment (YDSE) using either photons or 
electrons, quark flavour oscillations in the neutral kaon system and neutrino oscillations. 
In all cases the total probability amplitude for the system, Api, decomposes into the sum 
of two probability amplitudes, Ap I and Apf. 

A FI = A FI + Apj = {Apjle^ + \A^\e^ B (9.1) 

For the YDSE experiments, Apj and Apj describe experiments where either one of the 
slits is closed. For the neutral kaon or neutrino oscillation experiments, these amplitudes 
correspond to flavour mixing scenarios where only a single mass eigenstate is produced. 
The experimental observable, Pfi, is the probability to observe one of the set of states 
F: \f m ), m = 1,2,... given a prepared state in the set I: \ii), l = 1,2,... (see Section 2 
above). According to Feynman's principies I and III ((2.1) and (2.3)): 

Pfi = \A FI + Afj\ 2 

= \Apj\ ~\~ \Apj\ + 2Re[Apj Apj\ 

= \A FI \ 2 + \Af I \ 2 + 2\Af. I \\A* I \cas(<l> B -<j> A ) (9.2) 

The interesting physical phenomenon, in every case, is described by the interference term 
in (9.2) that depends on the phase difference <\>b—$a- The origin of the latter is different in 
all the four examples under discussion. In a YDSE, the phase difference is a consequence 
of different path lengths corresponding to the two probability amplitudes just as for 
the Michelson interferometer discussed in the previous section. However the physical 
mechanism that actually generates the phase difference is quite different in the two cases. 
For the photon experiment it is the space-time propagator of the excited source atom. For 
the electron experiment it is the space-time propagator of the electron itself. For neutral 
kaon oscillations the path difference is constant and the phase difference is a result of the 
different space-time propagators of the K5 and Kl mesons. The propagators are different 
because of the different masses of these states. As previously pointed out [11, 12, 13], 
for the case of neutrino oscillations, both the propagator of the unstable source particle 
whose decay produces the neutrino, and the propagators of the neutrino mass eigenstates 
Ui and v 2 give important contributions to the phase difference. In each case a realistic 
experiment is described in which all relevant physical parameters are discussed. In most 
previous discussions only simplifed models have been used. These models give essentially 
the same results as the path amplitude calculations for both the YDSE as well as for 
neutral kaon oscillations, but a markedly different one for neutrino oscillations. 
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Young Double Slit Experiment with Photons 




Figure 10: Geometry of a Young double slit experiment. Partides are produced at the 
source S and detected at D. For the photon experiment, S is a single excited atom. For 
electrons S is the exit of an electrostatic accelerator producing a narrow electron beam 
with the Gaussian momentum profile of (9.13). 



A schematic layout of a typical YDSE is shown in Fig.10, where the important geo- 
metrical parameters are defined. The photon source is a single excited atom at S. The 
photon passes through the slits A or B in correspondence with the probability amplitudes 
Ads^d) 1 or Apgito) 1 , and is detected at D at time to- Fig.10 shows the projections of 
the alternative paths of the photon in the xy plane, which is perpendicular to the long 
edges of the slits. In order for the photon to arrive at D it must be Fresnel diffracted, in 
the xy plane, at one of the slits. It is assumed that the height, h, of the slits is sufficiently 
small for this to occur, and also that the y-dependence of the detection rate, resulting 
from diffraction (i.e. due to integrating the contributions corresponding to paths passing 
at different positions within the slit) may be neglected for small vàlues of yo (see Fig. 
10). The widths of the slits, w, are assumed to be sufficiently long that the photon follows 
an essentially rectilinear trajectory in the xz and yz planes. The probability amplitudes 
^ds(^-d) 7 an d ^ds^dV are given directly by (4.8) above (compare Fig.l and Fig.10): 



Ai*(t D y 



l DS 



A 

t'ta 

À 



Adiffhwexp 



Adiffhw exp 



* < 


i 










i ( 


. 1 


n {KC 






~*2^ 



ta + t 



tb + r 



(9.3) 
(9.4) 



Neglecting the angular dependence of Adiff and setting = r% = l in the denominator of 
the RHS of (9.3) and (9.4) (a good approximation for small vàlues of yo), substituting the 
path amplitudes in (9.2) and integrating, with suitable limits, over to, as done above the 
the similar case of the Michelson interferometer, gives the result, similar to (8.7) above: 



P FI = P FI {t D < oo) = 2r s \Al\ 2 {l + exp 



(?b ~ r A ) 
2ct s 



cos k(t b - r A )} 



(9.5) 



where 



A 1 = 
a q — 



iAnhw 
2irlr' 



(9.6) 
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Assuming that yo, d <C l, it follows from the geometry of Fig.ll that: 



Ar = r B - r A 

giving in (9.5): 

P FI (y D ) = 2r s \Al\ 2 



2{d+\)y D 



(9.7) 



1 + exp 


(d+fW 


CT S l 





COS 



2 2 (2it{d+\)y D 
~ 4r 5 |^pcos 2 - 



(9.8) 



where in the last line of (9.8) the small exponential damping correction to the interference 
term has been neglected. This is identical with the well-known result derived using the 
classical wave theory. The spacing, Ay D , between adjacent interference fringes is give by 
(9.8) as: 



For typical experimental vàlues: l = lOcm, d + h/2 = lmm, À° = 5893À, Ay D is 29 
fiïn. For the Sodium-lines Tg at = 5.4ns, so that the damping correction in (9.5) is only 
exp(— 1.8 x 10~ n n) for the nth interference fringe. 

Young Double Slit Experiment with Electrons 

The actual realisation of Feynman's gedankenexperiment, discussed at length in chap- 
ter 1 of Vol III of [3] was not done until the early 1980s [50]. A closely analagous exper- 
iment was performed by Tonomura et al [49]. In this electron biprism experiment, the 
double slits of a YDSE are replaced by a region of elèctric field around a thin wire that 
provides an equivalent deflection of the electron trajectòries to that provided by Fresnel 
diffraction at the slits in a YDSE. 

In electron quantum interference experiments the electrons are typically produced by 
a 'direct current' source using thermionic of field emission. There is thus no source with 
a well-defined quantum state yielding a known time-dependent phase as in the optical 
experiments considered above, or the neutrino oscillation experiments to be discussed 
below. 

Referring to Fig.10 and using (3.11) for the space-time propagator of the electron, the 
probability amplitudes analagous to (9.3) and (9.4) for the electron case are, on making 
similar geometrical approximations to those made in (9.5) above: 



A A e 



aB e 
^-DS 



= *f 
= *if 



exp 
exp 



■-[E e At A - Pe (r' + r A )\ 
~[E e At B - Pe (r' + r B )} 



f(Pe)dp e 
f(Pe)dp e 



where 



_ A D [3 2 e A diff hwA s 

A °- Vi 



(9.10) 
(9.11) 

(9.12) 
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In these equations, As is the electron production amplitude, Adiff, the amplitude for 
difïraction at a slit and A D the amplitude for the electron detection process. The function 
f(p e ) is the normalised distribution of electron momenta at the detector. In the following 
a Gaussian form will be used: 

f(p) = —= — exp 



(9.13) 



It can be seen that the two important differences between a photon and an electron 
YDSE. are that, (i) in the electron case the phases of the probability amplitudes are 
determined by the space-time propagator of the electron, whereas in the photon case 
the phase is determined entirely by the coherent source, the contribution from the pho- 
ton propagator vanishing, and (ii) unlike on-shell photons, the electrons do not have a 
constant velocity, so that their actual velocity or momentum distribution must be taken 
into account in order to construct the probability amplitudes (9.10) and (9.11). Because 
of the photon's constant velocity it must be produced at different times in amplitudes 
corresponding to paths of different lengths. However, the electron can be produced at 
the same time in both probability amplitudes (i.e. Aí^ = At B in (9.10) and (9.11)) and 
produce the same detection event, at a well defined time, due to different velocities along 
each path. In this case, unlike for the photon YDSE above and neutrino oscillations, to 
be discussed below, there is no contribution to the interference phase from the production 
amplitude As- 

The Lorentz-invariant character of the phase of the space-time propagator of a particle 
(see (3.11) above) enables the phases of the complex exponentials in (9.10) and (9.11) to 
be written in the following four equivalent ways: 

mc 2 r mc 2 t (mc 2 ) 2 r (mc) 2 r . „ 

ó = = = -- — = -- — (9.14) 

r h -fh Evh P h y ! 

Where 7 is the usual relativistic parameter l/Vl — fi 2 -, fi = v/c = pc 2 /E = r/(ct). 
Substituting the appropriate vàlues of r and p then gives for the phase difference in (9.2): 

^_^ A = ^írl+iA.(il±isï) (9 . 15) 

c 2 h \ p A Pb ) 

For equal production times of the electron in the paths A and B the condition: 

+ + + r ' ' + rA r ' + r B / Q1R x 

t = t A = t B = = (9.16) 

va v b 

is respected. It is shown in Appendix F that, in this case, the RHS of (9.15) may be 
written, to first order in the small quantity Ar/r' = (r B — r a)/t' , as: 

pAr 2nAr 

<t>B ~ ^ = = -JDB ( 9 - 17 ) 

where 

__Pa+Pb , a . Q x 

p= (9.18) 
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and 

Xf B = - (9.19) 
P 

Which is de Broglie's formula for the quantum wavelength of a particle. As further 
discussed in Appendix F, the equal time condition (9.16) is essential to obtain this relation 
in the Feynman path amplitude approach. The relation (9.16) also implies the following 
velocity difference for the electron for the paths A and B: 

v B -v A = C -^p^( PB - P a) = ^ + 0[(Ar) 2 ] (9.20) 

where E,v and r are defined similarly to p in (9.18). 

Taking the modulus squared of (9.10) or (9.11), and performing the momentum inte- 
gration, gives: 



>A,B ei2 _ i Ae\2 
V DS — l^ol 



f J e i{<PA, B {p)-<t>A,B{P')) f(p)f( P ')5( P ~ P ')dpdp' 

Al\ 2 (9.21) 



2y / 7T(Tp 

The equal time condition gives, for the interference term in (9.2): 

Iab = 2\Aj I \\A^ I \cos((f) B - <p A ) 

= 2\A e \ 2 Re \ / / e ss f(p A )f(p B )5(p B ~Pa~ Ap)dp A dp B 



Re 



na 2 p 



/ e R e 2ct p e dp A \ (9.22) 



where, from (9.20), 



Ap = p B -p A ^p(^) -^- + 0[(Ar) 2 ] (9.23) 



The p A integral is readily evaluated by the method of 'completing the square' to yield the 
result: 

The first exponential damping term takes into account the fact that a non-vanishing 
spread in electron momentum is necessary for interference to occur - a consequence of the 
imposed condition of equal production times for the electron in the two paths. The second 
exponential describes damping due to the width of the momentum distribution. If it is 
very wide the probability to satisfy the condition (9.20) becomes very small. Inserting 
(9.21) and (9.24) in (9.2), as well as the expression (9.7) for Ar, gives finally: 



PfiÍVdY = 



\A 



e 12 

ol 



p 



Ap 2 a p Ar 2 ( 2(p)(d + h/2)y D \ 



l + exp{-[(-f) 2 + K— Y]}cos 



2(7/ v 2h ' JJ \ m J 



* 2 m m >( W + W»° ) (9.25) 
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where in the last line the typically small interference damping corrections 25 have been 
neglected. It can be seem by comparing (9.8) and (9.25) that, in this approximation, 
photons and electrons with \® B = A° are predicted to produce identical interference 
patterns in geometrically indentical YDSE. 

Strange Quark Flavour Oscillations 

In this case the Feynman path amplitude description of a specific experiment [51] will 
be considered. This is done to ensure that the effects of all possibly relevant physical 
parameters are properly taken into account in the discussion. In the experiment, Ks or 
K L mesons were produced via the processes: n~p — > A(K S ,K L ) by a 1.01 GeV/c pion 
beam. The initial state proton was either free within a chemically bound hydrogen atom, 
or bound in a carbon nucleus, of the polythene/plastic scintillator target. Roughly 50% of 
the interactions occured on bound protons. Semi-leptonic decays Ks, K L — > 7r ± e T z/ were 
observed, and their proper time intervals calculated. In this experiment the paths A and 
B in (9.2) correspond to the space-time propagation of the K s or K L mesons respectively, 
so in the following the labels A,B are replaced by S,L. It is assumed that either meson is 
detected at a fixed distance, L, from its production point. 

The detailed physics underlying the production amplitudes of the mesons is not well 
understood, but is not important for the description of the flavour oscillation phenomenon. 
An ss quark pair is produced by the strong interaction. The s quark is bound in the A, 
while the s undergoes a flavour- changing weak charged-current interaction which produces 
one of the mass eigenstates |K S ), |K L ) of the neutral kaon system. These states contain 
within their quark substructure both sd ( |K )) and sd ( |K )) components. Neglecting 
small CP-violating contributions, |K S ) and |K L ) are eigenstates of CP: 

|K S ) = -J=(|K ) - |K )) CP = +1 (9.26) 

|Kl> = ^(|K ) + |K ») CP = -1 (9.27) 

where the states |K ) and |K ) are related by the CP operator as: 

CP|K ) = -|K ), CP|K ) = -|K ) 

The s quark produced in association with the A couples only to the |K ) component of 
|K S ) and |K L ). It then follows from (9.26) and (9.27) that, neglecting CP-violating effects, 
the |Ks) and |Kl) production amplitudes are equal: 

(K s A|T|7r» = (K l A|T|tt» = A P (9.28) 

The decay processes Ks,Kl — > 7r ± e =F z/ have been found, experiment ally, to respect the 
'AS 1 = AQ Rule', which is predicted in the Standard Electroweak Model, by consideration 
of the W-boson exchange diagrams that mediate these transitions. This rule predicts that 
the |K ) components of |K S ) and |K L ), decay only into the channel 7r - e + z/ while the |K ) 

25 For example, the damping corrections for the nth interference fringe of the experiment of Refer- 
ence [49] where p = 229MeV/c, <r p /p = 6.0 x 1CP 7 (length of wavepacket 1/xm) and r' + r ~ 2m are: 
Ap/{2a p ) = Tvfhl[o v (r' + r)\ = ní.7 x and a p Ar/(2h) = mrcr p /p = nl.9 x 
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components decay only into n + e v. The equations (9.26) and (9.27) then lead to the 
following relations between the semileptonic decay amplitudes: 

(7i-e + u\T\K s ) = (7r-e + z/|T|K L ) = -(7i + e-u\T\K s ) = (ir + e~ v\T\K L ) = A eK (9.29) 

On the assumption that the Kg, Kl mesons are produced at space-time points xps, xpl 
respectively, and decay at xp, the probability amplitudes for the production of e + , e~ via 
propagation of Kg, Kl are given by (2.2) and (2.3) as: 

A e + S (K S ) = A eK J (D\K s \PS)f(ps)dp s (9.30) 

< 5 (K L ) = A eK I (D\K L \PS)f(p L )dp L (9.31) 

A e DS (K s ) = -A eK J (D\K s \PS)f(p s )dp s (9.32) 

A e DS (K L ) = A eK J (D\K L \PS)f(p L )dp L (9.33) 

The space-time propagators of the mesons are given by (3.12) as: 

[m l -tT i /(2c 2 )]c 2 Ar, 



(D\K i \Pt) = ^exp 



i = S,L (9.34) 



The amplitudes f(ps), ÍÍPl) describe the momentum distribution of the mesons. As 
discussed below the dominant source of momentum smearing is radiative corrections rather 
than variation of the physical masses Wl and Ws, of order T L /c 2 , T s /c 2 , of K L and Kg, 
due to their unstable nature. It is shown below that the precise form of f(p) does not 
affect the final result for the detection probability. For convenience, f{p) is normalised so 
that 

f(p) 2 dp=l (9.35) 



o 

Since the decay widths are relatively large, indeed T s = 2.1(m L — m s ), the kinematical 
effects of the variation of the physical masses of the mesons are of a similar size to those, 
generated by the difference of pole masses, that underly the whole 'flavour oscillation' 
phenomenon and so cannot be, prima facie, neglected. This variation of the physical 
mass, W, modifies the exact expression, (9.14), for the propagator phase in the following 
way: 

mc 2 T mc 2 t (mc 2 ) 2 L mc 2 WL 

^ = h~ = ~^h~ = ~E(w)v{w)h = ph~ (9 ' 36) 

where 

E(W) ee ^(Wc 2 ) 2 + (pc) 2 , v(W) = ^- ) 

The phase then depends, in general, on both the physical mass and the velocity (or 
momentum) of the particle. Unlike that of the photon in physical òptics the particle 
velocity is variable, and unlike the electron in the YDSE discussed above, the velocity is 
not fixed by the value of the momentum. 

In order to proceed further it is necessary to discuss the space-time structure of the 
production and detection events corresponding to the probability amplitudes. Only the 
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spatial distance L between the production and detection events is necessarily constant 
in the probability amplitudes (9.30)-(9.33). To give some feeling for the variation of the 
kinematical quantities and time intervals involved one may ask the following qüestions: 
(1) What is the difference of momentum necessary to compenstate the change in velocity 
due to the K S -K L mass difference? (2) In case of equal momenta for the K s and K L 
what is the difference in their production times in order to arrive simultaneously at the 
detection event? If the velocities of the Ks and Kl are equal, their momenta pl and ps 
must satisfy the condition: 



vs l v L 



c . i-\ + (±L·l\2 r . /i + (1L·l\2 

that is 



(9.37) 



(9.38) 



(9.39) 



W s _ ps 
W L ~ p L 

so that 

<M = Wl) - (W s ))c = (m L - m s )c = ^ ^ 1Q _ 14 
(p) (p) (p) 

where the value (p) = 194 MeV/c, corresponding to the centre-of-mass momentum of the 
Ks and Kl in the experiment [51], and — ms = 3.49 x 1CT 12 Mev/c 2 have been used. 
This is many orders of magnitude smaller than than the momentum smearing of the K s 
and K L due to initial state photon radiation. This is estimated, using the soft photon 
radiatve correction formalism of [52] to be: 

(Sprad) , w 1n _ 2 



<P> 



= 4.2 x IÓ' 2 (9.40) 



for the experiment [51]. In these circumstances the path amplitudes in (9.30)-(9.33) 
for different momenta and velocities, within the appropriate range given in (9.39), must 
contribute equally to the sum over intermediate states in (2.3), so that no damping of the 
interference term is to be expected from momentum or velocity variation. This is verified 
in the following calculation. 

To answer the question (2), the difference of production times: AtsL = tps — tph for 
Ks and Kl of fixed momentum, in order to arrive simultaneously at the typical distance 
cP'jTs from their production point is calculated: 

M SL = {mí ~™ s)Ts (941) 
bj 

Some vàlues of AtsL for different assumed vàlues of p and E are presented in Table 2. 
Associating these with the lifetime of a hypothetical unstable source particle gives decay 
widths of such a particle in the range from 1 to 200 MeV. whereas the 'characteristic time' 
26 of the strong/weak interaction process that produces the AKs or AKl systems must be 
much shorter than this as no such resonant state is actually produced. In the following 



26 This is an interesting concept that has not been addressed, to date, by any physical theory. Presum- 
ably it takes some non-vanishing time for the quark rearrangement and creation processes that convert, in 
the present example, the ir~p system into AKs or AKl, to occur. The possibility to observe such a time- 
interval in decay processes by measurable deviations from the exponential decay law or the Breit-Wigncr 
line shape of an unstable particle has been suggested [53]. 
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p (GeV/c) 


0.01 


0.1 


1.0 


10.0 


100.0 


E (GeV) 


0.498 


0.508 


1.117 


10.1 


100.0 


At SL (sec) 


6.27 x 10" 22 


6.14 x 10" 22 


2.8 x 10~ 22 


3.1 x IÓ" 23 


3.1 x 10" 24 



Table 2: Momentum and energy dependence of At S L, the difference in production times 
of K s and K L of equal momentum in order to arrive simultaneously at distance c(5^t s 
from their production point. 



it is assumed that the 'characteristic time' is so short in comparison with the vàlues of 
AtsL shown in Table 2, that the production time of the Ks and Kl in the probability 
amplitudes (9.30)-(9.33) is the same in every case, i.e. tps = tph = tp- Thus, since 
the detection event occurs some defmite time, the K s and K L are assumed to have equal 
velocities in the alternative probability amplitudes. Note here the difference with the 
electron YDSE. There also the equal production time assumption is made, leading to the 
conventional de Broglie wavelength for the associated 'matter wave', but in this case due 
to the unique physical mass of the electron, and the different path lengths in the YDSE, 
different velocities are necessary. 

Substituting (9.34) into (9.30)-(9.33) and using the expression for ATj given by the 
second and last members of (9.36), it follows that: 



tàKs, L )| 2 



I^ds(Ks,l)| 2 



A eK [3Ap 



exp 



f(p's,L)f(PS,L)S(p S ,L - (Ps,L)dp's,LdPS,L 

^s,LW s ,Lfn s ,LL 



\A 



ep | 



exp 



(p)h 
^s,LW s ^ L m SiL L 

{ P )h 



(9.42) 



where the approximations (3 L ~ f3 s = (3 and W L ~ W s — m = (m L + m s )/2 have been 
made and (p) denotes the average value of pl or p s . 

Replacing the path labels A and B in Eqn(9.2) by S and L respectively, the interference 
terms I e SL for the detection of e ± are given by (9.30)-(9.36) as : 



Isl ~ i2i?e \A ep \ 



II 



exp 




T s Ws | T L W L 



Ps Pl 
xf(Ps)f(PL)S(W s p L ~ W L ps)dp L dp s 



exp 



ic 
T 



m s W s m L W L 



Ps 



Pl 



(9.43) 



where the 5-function imposes the equal velocity condition (9.38) and the approximation 



2 

PsPl — P has been made. Performing the ps integral in (9.43) gives: 

; /exp 



Ig L — ±2i?e \A ep \ 



(T s + V L ) W L L \ f ic 2 Am LS W L " 
1 exp < : L 



2h 



Pl 



Pl 



Xf(^ L· )f(PL)d PL 



(9.44) 



where Am LS = m L — m s . Since the ratio Wl/Ws differs from unity only by quantities 
of order T s /c 2 m ~ (3.49 x 10- 12 MeV)/(498MeV) = 7.0 x 10~ 15 , then, to a very good 
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approximation, the replacements: Wl/Ws = 1, Wl = va may be made in (9.44) giving: 



1 SL 



±2Re \A ep \ exp 



(r s + r L )mL' 



' me 2 ArriLsL" 



2h (p) , 

where the normalisation condition (9.35) has been used. 



cos 



(9.45) 



Subsitituting (9.42) and (9.45) into (9.2) gives finally for the probability, P(e ± , L, (p)) 
to detect e ± at distance L from the production point of Ks, Kl with momentum (p): 



P(e ± ,L, (p>) 



ep] 



exp 



Tsmc 2 L 



±2 exp 



h( P ) 

(F s + T L ) mc 2 L 



+ exp 



2h 



cos 



h(p) 

'mc 2 AmLsL^ 



(9.46) 



On the assumption that the Ks and Kl mesons have the same velocity this probability 
may be expressed in terms of the proper time interval r between their production and 
decay via the relation (see (9.36): 



WL mL 



(p) ~~ (p) 



(9.47) 



so that 



P(e ± ,r) 



\A 



eP\ 



exp 



h 



+ exp 



±2 exp 



(Tg + T L )r 

2h 



cos 



rvr 

h 

^ Am LS c 2 T 



(9.48) 



This is just the formula that has been used till now to intèrpret such quark flavour 
oscillation experiments [51, 54, 55]. Thus the Feynman path amplitude calculation, where 
careful account is taken of all relevant physical parameters, particularly with respect to 
the space-time structure of the production and detection events, gives the same result as 
a simple calculation of the phase difference (ps — neglecting all such considerations, 
using the time- dependent Schròdinger equation in the respective rest frames of the Ks 
and Kl mesons.. 



Neutrino Oscillations 



Previous discussions of neutrino oscillations using the Feynman path amplitude for- 
malism can be found in [11, 12, 13]. The experiment to be discussed here is detection of 
the process: un — > e~p, where the incoming neutrino is produced by pion decay at rest: 
7r + — > The result found generalises in a straightforward manner to any neutrino 

oscillation experiment where the neutrino source is at rest. The experiment is shown in 
more detail in Fig.ll. A slow 7r + comes to rest at time ío m a stopping target T (Fig.ll 
a)). For simplicity, the case of only two neutrino flavours is considered. The neutrino 
mass eigenstates \ui) and \v 2 ) have pole masses m\ and m 2 where m\ > m 2 27 . Fig.ll 

27 Note that the subscripts '1' and '2' are arbitary, and do not necessarily correspond to fermion gen- 
eration number. In the presently favoured interpretation of atmospheric neutrino oscillations [58] , where 
there is an important contribution from pion decay neutrinos, the mass eigenstates are those associated 
with the second and third generations 
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b) and c) show two alternative histories for the stopped pion. In Figll b), the decay 
7r + — > [x + vi occurs at time íi, in Fig.12 c) the decay n + — > yU + z/ 2 occurs at a later time t 2 . 
With a suitable time difference t 2 — 1±, the neutrinos corresponding to the two alternative 
paths arrive at the detection event at the same time, tn and are thus indistinguishable 
(Fig.ll d)). The corresponding probability amplitudes A^, then add, as in (9.1), to 
give the total probability amplitude for the experiment. The labels 'A' and 'B' in (9.1) 
and (9.2) are here replaced by '(1)' and '(2)' corresponding to the mass eigenstates v\ and 

Note that the physical situation is strictly analagous to that of the Michelson Interfer- 
ometer or the photon YDSE. In these experiments the photon must be emitted at different 
times by the source, in the interfering amplitudes, in order to arrive in time coincidence 
at the detection event. In these two experiments the time difference is necessary to com- 
pensate the non-equal path lengths. In a neutrino oscillation experiment, both paths have 
the same length, and the time difference of neutrino emission is necessary to compensate 
for the different neutrino velocities. 

Introducing the process amplitudes for the 7r-decay and neutrino detection processes, 
and the space-time propagators of the pion source and the neutrinos according to (3.12), 
the probability amplitudes for propagation of the states \v\) and \v 2 ) are: 

A% = {e~p\T\un)U ei {D\u l \€)U l·xl {u l· i + \T\Ti + ){i\Ti + \Q) (i = 1,2) (9.49) 

the space-time propagators of the source pion (assumed to be at rest so that At = At) 
and the neutrinos are: 



(i|7T + |0) = exp 
(D\uí\í) = exp 



. {m n -iT n /(2c 2 )]c 2 (U-t ) 



(i = 1,2) (9.50) 



. mjC 2 (T D - Tj) 

h 



fi 

(i = 1,2) (9.51) 



The process amplitudes for pion decay and neutrino detection are written as products 
of 'reduced' process amplitudes and elements of the Maki-Nakagawa-Sakata [57] matrix 
Uu that describes the strength of the charged-current coupling between lepton i and the 
neutrino mass eigenstate |z/j): 

(e-p\T\u in ) = U ei (e'p\T\un) (i = 1,2) (9.52) 
{v ll· i + \T\u + ) = U^{u l· i + \T\u + ) (i = 1,2) (9.53) 

In the case of only two charged lepton and neutrino flavours, all the elements of the MNS 
matrix are real and are specified by a single mixing angle 612'- 

U e i = U e2 = cos0i2, U e2 = -C/ M i = sin 0i2 (9.54) 
Denoting the phase of as <p {i \ (9.49)-(9.51) give: 



c 2 



mi{r D - n) + m n (U - t Q ) (i = 1,2) (9.55) 
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C A T C b 



< /d > 



b) 




(71+) 




l 1 

► 










Vi 



c) 




(71+) 




t 2 


^2 








► 

V 2 



d) 




Figure 11: An experiment to measure neutríno oscillations following pion decay. a) a n + 
comes to rest in a stopping target T at time tç,. b) and c) show alternative histories of 
the stopped pion. In b) the decay 7r+ — > occurs at time t\. In c) the decay 7r+ — > /iz/ 2 
occurs at time t 2 - If mi > m 2 , and with a suitable value of í 2 — ti, the neutrinos may 
arrive at the same time to in the alternative histories at the the detection event, as shown 
in d). Since the detection event (Vi, z/ 2 )n — > pe" does not distinguish the neutríno fiavour, 
the amplitudes corresponding to the different histories in b) and c) must be added, as 
in (9.1), to give the probability amplitude for the experiment. Note the similarity of the 
sequence of space-time events to that in the Michelson interferometer of Fig.8. 
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Using the last member of (9.14) to re-write the first term on the RHS of (9.55) and noting 
that: 

U-t = t D -t - - (i = 1,2) (9.56) 



gives 



h m w (t D - t ) 

Pi Vi 



Po 



2po 



+ m n (t D - t - -) + 0[mt] (i = 1,2) (9.57) 

c 



where the neutrinos are assumed to be ultra- relativistic so that: 



Vi = c 1 - 



c 2 m 2 " 



2pg 



1,2) 



(9.58) 



In (9.57) and (9.58) p is the centre-of-mass momentum of a massless neutrino produced 
in the decay n — > /lis 



Po = ml -ml 

C 2m n 



= 29.8 MeV/c 



(9.59) 



The replacements A F 'j — > A^}^ in (9.2) give then, for the probability of electron detec- 
tion at distance L from the neutrino source: 



efj, 



roo 

/ dt D 

J tvrt in 



(e-p\T\isn)(isfi + \T\iT + ) 



x \ |C/ e i^i| 2 exp 



r /■> 

1 n L·i 



í'l 



h 



exp 



+ |^e2^"u2| 2 exp 



^(to — to) 



h 



r t /z 

% 



- c 2U el U e2 U l·l iU fi2 exp 



2^ 



cos 



'A(m 



12; 



,2^2 



cm,, 
2^o" 



- 1 



f) 1(9.60) 



where 



A(m 12 ) 2 = m l _ m \ 



In (9.60) the neutrino times-of-flight, t{\ have been introduced: 



t{ 1 = t D - U = - = - I 1 + 



c 2 m 2 ' 



2pg 



+ 0[mf] (i = 1,2) 



(9.61) 



and the unobserved detection time to has been integrated out. Consideration of Fig.12 
shows that the vàlues of í m j„ are: ío + íi for the propagation of neutrinos of mass mi only, 
to + í 2 f° r the propagation of neutrinos of mass m 2 only, and í + 1{ 1 for the interference 
term in (9.60). This is because both neutrino paths must be possible if interference is to 
occur, and t{ 1 > t{ 1 . Performing the tn integrations in (9.60) and substituting the vàlues 
of the MNS matrix elements from (9.54) then gives the final result: 



2h 

*en — -p 



(e p\T\un) (z//i + |T|7r + ) 



x < 1 — exp 



L 

IVA(m 12 ) 2 L 
Ahpl 



sin 2 9 12 cos 2 9±2 
A(m 12 ) 2 c 2 



cos 



Po 



cm w 
2p7 



- 1 



h 



(9.62) 
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The exponential damping correction of the interference term due to the non-vanishing 
pion lifetime is completely neligible 28 

The prediction (9.61) is readily generalised to the case of an arbitary two-flavour 
neutrino oscillation experiment with a stationary source,^: 



to yield the prediction: 



X s + U + Vi + T D ^ X D í k 



kj 



2h 



(£ k X D \T\uT D )(u£,X s \T\S) 



x 



{\u kl u n \ 2 + \u k2 u j2 \ 



+2U kl U k2 U jl U j2 cos 



A(m 12 ) 2 / R m 



m s 



1-RIJ h 



(9.63) 



Here li {li) denotes a charged lepton (antilepton) of fermion generation i and 

R m = {cm R ) 2 = {p{X s ) + pÇlj)) 2 

ms 

Here p{X s ) and p(Íj) are 4-vectors, m s is the mass of the source particle and m R is the 
effective mass of the particle, or particle system, X$ + lj, recoiling from the neutrino mass 
eigenstate. For the case of /3-decay, where: 

Ep = E{vi) + E(ïj) = constant 

the argument <f>^ — (f)^ of the cosine in (9.62) is replaced by [11, 12]: 



A{m 12 f 



Vu 



'El 
2Pu 



- 1 



(9.64) 



It is instructive to now compare in more detail the neutral kaon oscillation and neu- 
trino oscillation cases. In Table 2 is shown the difference in production times of Ks and 
Kl mesons of various fixed momenta in order that they arrive simultaneously, in the al- 
ternative histories corresponding to either Ks or K L propagation, at a detection event at 
a typical decay distance cf3^Ts- A similar comparison is now made for neutrino oscilla- 
tions by considering the difference of production time Aí i2 = t 2 — t\ in order that the 
neutrinos v\ and v 2 arrive simultaneously at the detection event, as shown in Fig.ll d). 
The distance, L{tt), between the production and detection events is chosen so that the 
phase of the cosine term in (9.62) is n rad. Thus: 

hpo 13. 8m 

W = A(m 12 ) 2 c 2 (^ -2) = A(m 12 ) 2 c 4 (eV) 2 ) (9 ' 65) 



28 Setting A{mi2) 2 c 2 L/{pQh) to unity, so that the argument of the cosine in (9.62) is of the same order, 
as necessary for an observable neutrino oscillation cffect, gives for the damping term: exp[— r,r/(4po)] = 
cxp[-4.0 x 10" 16 ]. 
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The corresponding production time difference is: 



where (9.57) has been used. Combining (9.65) and (9.66) gives: 

Aíai = 77^7 r = 8.22 x l(T 24 sec (9.67) 

Thus AÍ21 depends only on the kinematics of the neutrino production process and is 
independent of the neutrino mass difference. For neutrinos produced in pion decay, At 2 i 
is a factor of 3 x 1CT 16 smaller than the pion mean lifetime of 2.6 x 1CT 8 sec. This 
implies strictly equal amplitudes for the decay processes shown in Fig.ll b) and c) with 
completely negligible lifetime damping effects as discussed above. Further estimates of 
the latter may be found in [12]. The situation is thus dramatically different from the 
neutral kaon flavour oscillation case where, instead, the time differences for simultaneous 
arrival of equal momentum K s and K L mesons at the detection event shown in Table 2 
are much larger than the characteristic time of the production process. Thus the K s and 
K L are produced at essentially the same time. Because of the tiny K S -K L mass difference: 



Am 12 /m = 3.49 x 1(T 12 MeV/c 2 /497.7 Me V/c 2 = 7.0 



x 10 



,-15 



momentum and velocity smearing effects due to variation of the physical masses of the Ks 
and Kl and radiative corrections ensure that the probability amplitudes for Ks and Kl 
of the same velocity are closely equal, leading to the conventionally used formula (9.47). 
In the neutrino case however inspection of equal velocity condition analgous to (9.38): 

^ = * (9.68) 
m 2 p 2 

shows that it may be physically impossible for v\ and v 2 to have the same velocity in the 
interfering probability amplitudes. 

The current analysis of atmospheric neutrino oscillations where pion decays give an 
important contribution estimates the mass difference to be [58]: 

A(m 12 )V ~2x 10" 3 (eV) 2 (9.69) 

In order to respect the equal velocity condition (9.68) the smearing, Sp, of the neutrino 
momentum must be such that: 

5p _ Am 12 _ m 1 - m 2 ^ ^ 

p m rïï 

Now, 

A(mi 2 ) 2 = (mi — m 2 )(mi + m 2 ) = 2Ami 2 m (9-71) 
So, combining (9.70) and (9.71), equal velocities require that: 

5_p _ A(m 12 ) 2 

p ~ 2(m)2 [yJZ) 
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Assuming m = 0.5 eV/c 2 , which is consistent with all direct upper limits on neutrino 
masses [56], and inserting the experimental value of A(mi 2 ) 2 from (9.69) gives: 

| = 4x 1(T 3 (9.73) 
V 

The dominant source of momentum smearing in pion decay is radiation of real photons. 
Using the soft photon formalism of [52] to calculate the mean momentum smearing, (Sp ra d) 
due to soft photon emission gives: 

{Sprad) q r 1n -4 



<p> 



3.5 x 10" 4 (9.74) 



which is incompatible with the equal velocity condition (9.73). Indeed the latter condition 
can be consistent with (9.74) only if m is close to the current direct upper limit on m Ul of 
3 eV/c 2 [56]. In this case 5p/p ~ 1.0 x 10~ 4 , which is just compatible with (9.74). There 
is still, however, no physical reason why probability amplitudes with equal velocities and 
production times should be favoured over those with different velocities and different 
production times. 

The difficulty to justify the equal velocity hypothesis for neutrinos in view of the 
condition (9.68) has been previously pointed out in the literature [59, 60] but the authors 
of these papers still concluded that the 'standard' formula for the neutrino oscillation 
phase [58], which will now be discussed, was correct. 

As pointed out in [11], the derivation of the standard formula requires that the neu- 
trino velocities defined kinematically according to the formula Vki n = pc 2 /E are different 
for different mass eigenstates, but that the space-time velocities: v s -t = L/t are the same. 
These manifestly contradictory hypotheses (if neutrinos are indeed real (on-shell) partic- 
ules propagating in space-time) give, instead of the prediction of (9.63) for the oscillation 
phase: 

fe ^<,_^^(_^y i + K ] (9.75) 



the standard formula: 

A(m 12 ) 2 c 2 



0í w = _ 0(2) = L + 0[ f] (9 76) 

2p h 

This is actually the same result as found here for the neutral kaon oscillation phase (9.45) 
since: 

_ A ( mi +m 2 ) A(m 2 2) 
mAm 12 = (mi - m 2 ) = — - — - 



A simple derivation of (9.76) can be found in [58]. Each neutrino mass eigenstate is as- 
sumed to evolve temporally in its rest frame according to the time-dependent Schròdinger 
Equation. The corresponding phase increment is correctly given by the first member of 
(9.14) above. The Lorentz invariance of the phase is then used to write it in terms of 
laboratory-frame quantities: 

«^) .^^i!^) (i = 1>2) (9 . 77) 
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It can be seen that it is just at this point in the derivation, that the assumption of equal 
propagation times, t, for each mass eigenstate implies equal space-time velocities L/t 
whereas the kinematical velocities PíC 2 /Eí are assumed to be different. As discussed in 
detail in [11] it is this unphysical assumption in the approximation made in the last 
member of (9.77) that leads to the factor of 2 difference between (9.76) and the neutrino 
propagation phase given by the exact relation (9.14). The Standard formula (9.76) is 
derived from <f>^ — <j>( 2 \ as given by (9.77), on making the further assumptions that the 
neutrinos are ultra-relativistic, that L/t = c, and retaining only the leading A (777.12) 2 
neutrino mass terms. Evidently, since (9.77) assumes equal production times for different 
mass eigenstates, there is no contribution to the oscillation phase from the source particle 
in the Standard formula. 

The Standard formula (9.76) is still universally employed by the neutrino physics 
community. Claims that the Feynman path amplitude formula (9.75) is incorrect have 
been made [61] and rebutted [11, 62]. A critical discussion of the ad hoc 'wavepacket' 
approach, that is frequently used in discussions in the literature of the quantum mechanics 
of neutrino oscillations, that results, - from the point-of-view of Feynman's formulation 
of quantum mechanics embodied in the laws I-V presented in Section 2 - an artificial and 
unphysical blurring of the space-time structure of production and detection events, can 
be found in [11, 12]. 

The formula for the production time difference, (9.66) becomes, in the general case 
described by (9.67): 

Thus, when the mass of the system recoiling against the neutrino becomes very small in 
comparison with the mass of the source particle, At 2 i becomes very large, as does also 
the 'oscillation length, L osc , defined by writing the argument of the cosine interference 
term in (9.63) as 2ttL/L osc so that: 

This dependence of the oscillation length on the kinematics of the neutrino production 
process gives a straightforward experimental method to discriminate between the path 
amplitude prediction for the oscillation length, (9.79) and that predicted by the Standard 
formula (9.76) which is: 

L s í7 d = a ? P °i 5 (9-80) 



As proposed in [13], if a terrestrial long linebase 'z/^ disappearence' experiment such 
as K2K [63] is performed using neutrinos produced in the process: K — > fiu, rather 
than 7r — > fiv as in the existing experiment, the path amplitude formula predicts that 
the oscillation length is a factor of ~ 28 times longer for a kaon source than for a pion 
one yielding the same neutrino momentum. This implies a strong suppression of the í v il 
disappearence' phenomenon when using a kaon rather than a pion source. The Standard 
formula predicts the same oscillation length for pion and kaon sources giving neutrinos of 
the same momentum. 
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Experiment 



Photon YDSE 



Electron YDSE 



K Osc. 



v Osc. 



Sr = 0? 
St = 0? 
5v = 0? 



5<p s 



0? 



No 
No 
Yes 
No 
Yes 

-pAr 
1 



No 
Yes 
No 
Yes 
No 

-pAr 
1 



Yes 
Yes 
Yes 
Yes 
No 

A(rn LS ) 2 c 2 

2 (-2- 

\msc 



Yes 
No 
No 
No 
No 



A(m 12 ) 2 c 
m s 

2 (i 

msc V 



L 



-R 2 



Rn 



Table 3: Comparíson ofthe four 'two probabilty amplitude' experiments. The characteris- 
tics of each experiment are specifíed by the quantities: 5r, 5t, Sv, S<f) source and 5<f) part which 
are the diiïrences of; path length, time-of-üight, velocity, source phase and particle phase 
for the two interfering amplitudes. Also shown is the interference phase <pB — 4>a m (9-2) 
and the ratio of the effective wavelength defined in (9.80) to the de Broglie wavelength. 



The four different 'two probability amplitude' experiments that have just been dis- 
cussed are summarised in Table 3. Each experiment is first classified in terms of its 
space-time properties. These are specified by the differences between the path lengths 
(Sr), the times-of-flight (St) and the particle velocities (Sv) for the two paths. This in- 
formation is given as a Yes/No answer to the qüestions: Sr = 0?, St = 0? and Sv = 0?. 
Similarly the contributions to the interference phase from the source and the propagating 
partides are specified by the differences <50 source an d S4P art for the two paths and given 
as answers to the qüestions S<p source = 0? and 5<p part = 0?. In the last two rows the vàlues 
of the interference phase and the ratio /X DB are reported. The effective wavelength, 
\ e ff is defined for neutral kaon abd neutrino oscillations in terms of the phase difference 
Aip for paths of length L by the relation: 

- w\ (9 ' 81) 

while \ DB is the conventional de Broglie wavelength of the propagating particle as given 
by (9.19). 

The puzzling and difncult-to-understand concept of 'wave-particle duality' does not 
survive the detailed analysis of the superficially analogous, but in fact very different, 
experiments presented in Table 3. The 'wave' associated with the particle is different 
in every case. Although both the YDSEs may be simply (and correctly) interpreted in 
terms of the de Broglie wavelengths of the photon or the electron, the physical origin 
of the 'wave' is quite different in the two cases. For the photon it is the phase of the 
propagator of the excited source atom; for the electron it is the phase of the propagator of 
the latter. The effective wavelength for neutral kaons, although originating entirely from 
their propagators is quite different from that of the electron in the YDSE experiment. 
This is a result of the different space-time properties of the experiments as detailed in 
Table 3. For neutrino oscillations is again different and depends on the propagators 
of both the source particle and the neutrinos. 

The conclusion to be drawn from Table 3 is clear. In a strict, ontological, sense only 
the source and propagating partides exist as physical entities, and it is their particular 
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properties together with the space-time structure of the production and detection events 
that determine (ps — <Pa- The associated 'wave' is only a (sometimes useful, sometimes 
not) mathematical abstraction, as clearly stated long ago by E.J. Williams [64] 

The electron is, of course, a particle. The wave is in the mathematics. 

Although its inventor doubted it until the end of his days 29 the photon is also a parti- 
cle. The realisation of this is enough remove Einstein's mysterious and counter-intuitive 
postulate concerning the constancy of the velocity of light from the foundations of special 
relativity [66], and, by comparing photon properties with those of classical electromagnètic 
waves, to understand, in a new way, the bàsic concepts of quantum mechanics [67]. Also, 
all the mechanical aspects of classical electromagnetism may be quantitat i vely explained, 
at the quantum level, as an effect of the exchange of space-like virtual photons [68]. 



10 Summary and Outlook 



In this paper Feynman's path amplitude formulation of quantum mechanics has been 
applied to several specific problems. In particular many of the optical experiments dis- 
cussed in a qualitative way in Feynman's book 'QED' [1] have been calculated in detail. 
The problems of physical òptics addressed include: the calculation, from first principies, 
of the refractive index of a transparent medium, the laws of refraction and reflection of 
light at the interface of two transparent media, the rectilinear propagation of light in a 
homogeneous transparent medium, Fermat's Principle, the reflection coefficient of light 
at normal incidence from the interface between two transparent media, and temporal 
dependence of the fringe visibility in a Michelson interferometer. The path amplitude 
method is also used to analyse quantum interference effect s in Young double slit experi- 
ments (YDSE) using photons or electrons, in neutral-kaon quark flavour oscillations and 
in neutrino flavour oscillations. In all cases, except neutrino oscillations, agreement is 
found with the predictions of the classical wave theory of light or 'wave particle duality' 
(i.e. introducing the de Broglie wavelength of a massive particle in analogy with the 
wavelength of the photon), for the YDSEs, or, in the case of neutral kaon or neutrino 
flavour oscillations, solving the time- dependent Schròdinger Equation in the particle rest 
frame. 

All of the calculations presented are straightforward applications of Feynman's five 
laws of quantum dynamics, that provide a description of causally-related elementary pro- 
cesses that succeed each other in space-time. These laws, presented and discussed in 
Section 2, are: 

(I) The probabilistic interpretation of the Probability Amplitude (2.1). 

(II) The law of Sequential Factorisation (2.2). 

(III) The law of Quantum Mechanical Superposition (2.3). 
29 'What are light quanta? Today every scoundrel believes he knows the answer, but he is wrong.' [65] 
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(IV) The law of Composite Factorisation (2.4) 
(V) The Feynman Path Integral (2.9) 

The dynamics of the theory is entirely contained in V. The laws II- IV describe how 
the probabilty amplitude is constructed, while I gives its physical interpretation. The 
probability amplitude is built, according to II, as the product of a series of amplitudes 
describing elementary physical processes ordered in time. For the calculation of space- 
time dependent interference effects, the most important of these process amplitudes are 
the Green functions, of space-time propagators, of partides in motion or at rest. The 
derivation of these propagators from, respectively, the Feynman Path Integral and the 
time-dependent Schròdinger Equation, is described in Section 3. 

Before addressing specific problems of quantum interference, the physical interpre- 
tation of Feynman's formulation of quantum mechanics is discussed in some detail in 
Section 2. This interpretation is indeed quite different from that of previous formula- 
tions since the primary physical concept, the path amplitude, is neither a Hilbert-space 
state vector nor a wavefunction, and the concept of a particle, localised in space-time, 
as in classical physics, is also essential. Indeed consideration of classical propagation of 
partides localised in space-time is necessary for the correct calculation of the phases of 
the path amplitudes that are the physical basis of all space-time quantum interference 
phenomena. In this theory, partides, once created, propagate in space-time in a classical 
manner within each path amplitude. In the limit h — > the sum over paths is replaced 
by a single path, the one predicted by Hamilton's Principle of classical mechanics, and 
Feynman's laws I and II become a statement of Bayes theorem for the combination of 
conditional propabilities. 

Shortcomings of the conventional text book 'wave packet' representation of partides 
and the associated misinterpretation of the space-time Heisenberg Uncertainty Relation 
are pointed out. The relation of Feynman's formulation to other 'interpretations' of quan- 
tum mechanics: Consistent Histories, The 'Many Worlds' interpretation and de Broglie's 
Pilot Wave theory are discussed. The much debated 'measurement problem' of the Heisen- 
berg and Schròdinger formulations of quantum mecànics, seems, at least superficially, to 
be absent from Feynman's formulation. Perhaps new 'problems' remain to be discov- 
ered; but in any case, to the writer's best knowledge, no discussion of the 'philosophy' of 
Feyman's formulation yet exists. 

In Section 4 it is demonstrated that Feynman's path amplitude formalism, in which 
the dynamics is entirely contained within the Green functions of partides, predicts, in a 
completely general manner, the classical wave theory of light, under a certain well-defined 
condition. The latter is that the lifetime of the photon source must be much longer than 
the photon path differences, in the experiment under consideration, divided by the velocity 
of light. This is the case for essentially all quantum òptics experiments using a làser as 
light-source. All such experiments can then be correctly analysed using the classical wave 
theory of light, i.e. as done in well-known text books on òptics [33, 37], in terms only 
of spatial intervals and photon wavelengths, without taking into account the times of 
the events corresponding to different observed physical processes. Although Feynman's 
formulation predicts the existence of an effective 'wave field' U that satisifies Hemholtz' 
Equation (4.13), from which all predictions for interference and diffraction phenomena of 
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the classical wave theory may be derived, the physical source of the phase of U is solely 
the temporal propagator of the excited source atom, not the propagator of the photon 
itself. Thus any 'wave particle duality' in physical òptics is between the photon ('particle') 
and its source ('wave') not between the 'particle' and 'wave' aspects of the photon itself. 
Indeed the latter does not exist in the path amplitude formulation. Thus, starting only 
from the quantum properties of partides (their Green functions) their 'wave' properties 
may be derived in full. It can seen that the latter are then, albeit useful, mathematical 
abstractions that do not exist in the same strictly ontological sense as the partides. Given 
the a priori existence of partides and Feynman's rules for their quantum behaviour, 
all their so-called 'wave' properties may be derived. The inverse procedure is evidently 
not possible. The fundamental physical entities are then the partides, not the 'waves'. 
In consequence there is in fact no 'wave particle duality' in the physical description, 
though by chance, in some cases such as in Davisson and Germers's original 'matter 
wave' experiment [69] or the electron Young double slit experiment, a naive application 
of this idea gives the same prediction as the full path amplitude calculation. 

Finally in Section 4 it is demonstrated that the amplitude for forward diffractive 
scattering derived from the Helmholtz Equation by means of Green's theorem, can also 
be derived by direct summation of path amplitudes. The integration procedure used is 
similar to that employed in Section 5 to calculate the refractive index of a transparent 
medium in terms of the scattering amplitude of light by the àtoms of the medium. 

Section 5 contains three different calculations of the refractive index. In the first two 
it is assumed that the photon which interacts with the refractive medium is produced 
sufficently long after after the production time of the excited source atom that all ge- 
ometrically allowed paths of the photon are possible. The first calculation assumes a 
sheet of refractive material so thin that, to a good approximation, the refractive index 
may be calculated by considering a single scattering of the photon from the àtoms of the 
material. For this condition to be satisfied the sheet must be very thin -a small fraction 
of the photon wavelength. The second calculation removes this restriction and considers 
a block of uniform refractive material of arbitary thickness. In this case any number of 
scatterings of the photon from the àtoms of the material are taken into account. The 
formula obtained for the refractive index is the same (5.9) as that obtained in the single 
scattering case. The third calculation considers photons produced and detected promptly 
after the production of the excited atom. In this case, causality forbids photon trajectò- 
ries that deviate significantly from the classical straight line path between the source and 
the detector, reducing the number of possible interactions with the àtoms of the refrac- 
tive medium and therefore strongly reducing the refractive index -the so-called 'refraction 
annulment' effect. A simple experiment to detect this phenomenon (if indeed it exists) is 
shown in Fig.4. 

In Section 6 the laws of refraction, reflection and rectilinear propagation of light are 
derived by imposing the condition that the phases of the probabilty amplitudes of suitably 
defined experiments are stationary for variations of the spatial location of the photon path 
from its source to the detector. The generalisation of the stationary phase condition to 
the motion of an arbitary physical object gives the quantum mechanical explanation of 
Hamilton's Principle of classical mechanics. Since the phase of the probability amplitude 
is proportional to the effective propagation time of a photon from its source to the detector 
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Fermat's principle of least time is also derived. By considering the second derivatives of 
the phase of the probability amplitude, with respect to the positions and angles of the 
corresponding paths, the the size of the spread of photon trajectòries around the classical 
path corresponding to a stationary phase, is estimated. 

In Section 7 the experiment used in [1] to exemplify the application of quantum 
mechanics to òptics, the reflection of light at normal incidence from the boundary of two 
transparent media with different refractive indices is analysed in detail for the case of a 
realistic experimental set-up. It is essentially the same as the first quantitative quantum 
mechanical experiment ever performed -Newton's study of his 'Rings' [2]. The result 
for the reflection coefficient at a vacuum/medium interface, (7.8) is significantly different 
from the formula (7.10), first derived by Fresnel, that is to be found in all text books 
on òptics. For a refractive index n = 1.5, the reflection coefficient predicted by the path 
amplitude calculation is 0.028, some 40% smaller than the prediction, 0.040, of the Fresnel 
formula. Experimental discrimination between these two possiblities should be simple and 
straightforward. The present writer has been unable to find either in text books or in 
the published research literature any account of experimental verification of the Fresnel 
formula (7.10) for photons in the optical region of wavelength. 

Fresnel originally derived his formula for the reflection coefficient from an elastic-solid 
model of light waves. The same formula is obtained in classical electromagnetism by 
applying surface boundary conditions to the elèctric and magnètic fields of the different 
electromagnètic waves in refractive media of different refractive indices. The problem 
with this approach may be that although a classical electromagnètic wave is, certainly, 
a high density beam of monochromatic photons, a single photon is, equally certainly, 
not a classical electromagnètic wave, so it makes no sense to associate with it classical 
elèctric and magnètic fields. If the Feynman space-time picture of reflection is correct, 
the essential physics of the process relates not to the interface between media but resides 
in the amplitude A sca t for elàstic scattering of photons by the àtoms in the interior of 
the medium, as is manifest in the derivation of the formula (5.7) for the refractive index. 
This formula, and the formula (7.8) for the reflection coefficient are derived in exactly 
the same way by by integrating the path amplitudes corresponding to the scattering of 
the photon over all of the àtoms of the refractive material probed by the photon in the 
experiment. The interface occurs only geometrically as a limit on this spatial integration. 

This being said, it is certainly possible to envisage an experiment where the Fresnel 
formula would be appropriate and correct. For example, the reflection of a beam of 
microwave photons of sufficiently high spatial density from the surface of a dielectric. In 
this case the description of the beam as an electromagnètic wave will be appropriate, 
and the molècules of different dielectrics will be differently polarised, over macroscopic 
distances, by the elèctric fields associated with the microwave beam. The waves in different 
media will have different characteristics related by the boundary conditions on the elèctric 
fields of the waves, so that the Fresnel formula for the reflection coefficient should be 
correct, and the scattering amplitude of individual microwave photons from individual 
molècules less directly relevant. On the other hand, the macroscopic dielectric polarisation 
effects for a single optical photon entering a block of glass are clearly negligible and a 
purely quantum mechanical description is mandatory. For a general discussion of the 
different circumstances in which either the classical electromagnètic wave (high photon 
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density) or the quantum mechanical (low photon density) description is appropriate, see 
[67] 

The experiment described in Section 8, measurement of the temporal dependence of 
the fringe visiblity in a Michelson interferometer is one in which path length differences are 
deliberately chosen to be of the order of, or greater than, cts where Tg is the mean lifetime 
of the excited atòmic state that produces the observed photon. The classical wave theory 
of light, in which only spatial intervals are considered is therefore unable, in principle, to 
describe such an experiment. Thus new predictions specific to the path amplitude method 
are obtained. These predictions for the fringe visiblity as a function ot the time elapsed 
after the production of the excited source atom are shown, for different path length 
differences, in Fig.9. The path amplitude calculation of the effects of random source 
motion on the time-integrated fringe visibility found in Appendix E predicts that there 
should be no effect on the fringe visiblity from this source. This is in contradiction with 
the prediction, (8.14), of Rayleigh, based on the classical wave theory, that such 'Doppler 
broadening' effects should be large 30 . In order to test the path amplitude prediction of a 
vanishing effect due to source motion, experiments are required that clearly discriminate 
this effect from that of pressure broadening. So far this has not been done. 

In Section 9 analagous 'two probability amplitude' experiments in photon or electron 
òptics, neutral kaon oscillations and neutrino oscillations are analysed and compared. 
The understanding of the different results obtained, summarised in Table 3, requires 
careful consideration of the nature of the particle source and the partides themselves as 
well as the space-time structure of the production and detection events of the different 
path amplitudes. Only for the photon and electron YDSE does the simple 'wave particle 
duality' concept, that is the introduction of the de Broglie wavelength of the photon or the 
electron, give the same result as the path amplitude analysis. For the neutral kaon and 
neutrino oscillation cases, quite different effective wavelengths are found for the associated 
'matter waves'. For optical interference experiments the source of the 'wave' is solely the 
propagator of the excited source atom. As neutrino oscillation experiments (unlike the 
electron YDSE and neutral kaon oscillation experiments) have a similar excited source 
particle, it is clear, if the path amplitude description is indeed correct, that the source 
particle must give an important contribution to the oscillation phase also in this case. In 
addition there is the contribution of the propagators of the neutrinos themselves, which 
is analogous to the sole source of the interference phase in the electron YDSE and neutral 
kaon oscillation experiments. Before the work reported in [11, 12, 13] only this last 
contribution to the neutrino oscillation phase was taken into account. 

Although all optical and massive particle diffraction and interference experiments that 
have been performed to date and found to be agreement with the predictions of the clas- 
sical wave theory of light, or 'wave particle duality' can be interpreted, following the 
arguments given in Section 4 and 9, as evidence for the correctness of Feynman's formu- 
lation of quantum mechanics, the present paper contains several new predictions. If any 
of these predictions is not confirmed by experiment, Feynman's space-time formulation of 
quantum mechanics is not of general validity and answers must be sought to the question 

30 Indeed, as shown by the entries in the first and last rows of Table 1, the damping predicted by the 
Rayleigh formula is too large to explain Michelson's fringe visiblity measurements for the H r 3p-2s and 
Na D 3p-3s transitions. 
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why, in some cases, cited above, the predictions agree with experiment, and in others not. 
These new predictions are: 

(i) The refraction annulment effect (Section 5). 

(ii) The reflection coeficient of light at normal incidence (Section 7). 

(iii) Time dependence of fringe visiblity in the Michelson interferometer (Section 8) 

(iv) Absence of source motion damping of fringe visiblity in the Michelson interferometer 
(Section 8). 

(v) Different oscillation lengths for 'z/ M disappearence' using n — > jiu or K — > jiu as 
neutrino sources. 

It is clearly of interest to apply the methods developed in the present paper to other 
related physical problems. Some examples are: 

(A) The coeficients of reflection and refraction of light at non-normal incidence. For 
this it is necessary to take into account photon polarisation. This requires, for 
the path amplitude method, a complete quantum description of polarisation effects 
in both the photon production process (elèctric dipole radiation from a polarised 
source atom) as well as photon scattering processes (Rayleigh scattering of polarised 
photons). Another interesting related problem is the path amplitude description of 
bi-refringence. 

(B) Heavy quark flavour oscillations for entangled systems. Examples are: 

0^K L K s T(4S)^B H B L 

The latter is of considerable current interest due to the on-going experimental pro- 
grams of the BABAR and BELLE b-factories. Unlike the 'unentangled' process 
7i~p — > A(Kl, Ks) described above, both the <fi and the T(4S*) are coherent sources, 
similar to an excited atom in physical òptics. The corresponding vàlues of Aísl, 
At LH calculated using (9.41) are 6.12 x 10~ 25 s, 9.55 x 10~ 26 s repectively, where the 
value rriH — rn L = 3.28 x 1CT 10 Mev/c 2 [56] has been used for the neutral b-meson 
mass difference. These time differences may be compared with the mean lifetimes 
of the and the T(4S I ) source partides of 1.5 x 10~ 22 s and 4.7 x 10 _23 s [56], which 
are factors of 24.5 and 492 times greater than Aísl an d At LH . Thus, unlike for the 
case of n~p — > A(K L ,K S ) neutral K- or b-mesons can, in principle, have the same 
detection time in alternative histories with different mass eigenstates, provi ded that, 
as in physical òptics or neutrino oscillations, the partides are produced at different 
times in the alternative histories. In view of more stringent space-time restrictions 
imposed by the requirement of observation of both final state partides it is possible 
that the off-shell nature of the propagating partides is more important, in this case, 
in enabling interference, rather than production time differences or velocity smear- 
ing due to radiative effects. Only an actual calculation can show if this is the case 
or not. 
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(C) Effective refractive indices for neutrinos. It is important, in view of the great current 
interest in neutrino oscillation experiments to repeat the refractive index calculations 
of Section 5 above for the case of neutrinos interacting with matter. The model 
used till now to derive the effective refractive indices of neutrinos in bulk matter 
[70, 71] was constructed by analogy with that used to used to describe coherent 
regeneration [72] of neutral kaons for which (see Table 3 above) simple 'wave particle 
duality' is a good approximation to path amplitude calculations. However, in the 
neutrino case, the path amplitude calculation suggests that this analogy breaks 
down, and an important contribution to the interference phase originates in the 
propagator of the source particle. This effect has not, to date, been included in the 
calculation of neutrino refractive indices. 

One of the precepts of the Standard 'Copenhagen Interpretation' of quantum mechanics 
is that our direct exerience of the world can only be understood, in a rational manner, by 
an appeal to classical concepts 

The quantum theory is characterised by the acknowledgement of a 
fundamental limitation in the classical physical ideas when applied to 
atòmic phenomena. The situation then created is of a peculiar nature, 
since our interpretation of the experimental material rests essentially 
upon the classical concepts [9] . 

Is it indeed true that '...our interpretation of the experimental material rests essentially 
upon the classical concepts'? In the case of one 'experiment' that is performed by every 
human being 31 every day of her (or his) life without she (or he) being consciously aware of 
it, this statement seems to be quite untrue. The 'experiment' in question is our perception 
of colour. The photon detector D, forming an essential part of all the optical experiments 
discussed in Sections 2-8 above, was never precisely specified, but could well be one of 
the three types of cone receptor celis in the retina of the human eye. The detection 
efficiency of these receptors for photons is similar to that of a photo-diode with a flat 
spectral response in the optical region preceded by wide-band wavelength filters with 
roughly Gaussian acceptance profiles centered at À 7 ~ 440, 540 and 570 nm with full 
widths at half-maximum ~ 200, 200 and 300nm [73]. What we perceive as the colour 
'green' corresponds not always to some discrete or narrow range of vàlues of À 7 but 
in general to a weighted average (by the human nervous system) of the intensities of 
photons in different spectral regions that are incident on the three types of receptors. 
The 'interpretation of the experimental material' mentioned above by Bohr when this 
material is, say, the perception of the colour green or the colour red, varies enormously 
according to the circumstances, but in all cases this interpretation is possible only in 
terms of quantum, not classical, concepts. The sources of the photons for all the colour 
perception experiments that will now be discussed are excited àtoms in the Sun. Myriads 
of different spectral lines produced by spontaneous decay of these àtoms are responsible 
for the 'white light' arriving at the surface of the Earth in daytime. The experiments 
differ in exactly how certain of these photons come to impinge on the cone receptors of 
the retina. 

I step out on to my balcony and look at a geranium. I see that its leaves are green and 
31 Colour- blind people excepted 



72 



its flowers red. Why? The explanation is the different atòmic structures of the matter 
of the leaves and flowers. Photons of some wavelengths are absorbed, others re-emitted 
as fluorescent radiation by these àtoms. The re-emitted radiation, incident on the retinal 
receptors and suitably re-weighted gives the perception of green for the leaves and red for 
the flowers. In this case the physics explanation of colour vision is simple, the different 
parts of the geranium filter, in a different way, the sunlight photons that are incident on 
it. 

Nearby there is a shower of rain, and the Sun is lowering towards the western horizon. 
I look at the eastern sky and, at certain angles, see the green and red components of a 
rainbow. In this case the physics explanation of my perception of these colours is quite 
different. Now the photons responsible for the sensations of 'green' and 'red' have in each 
case closely similar wavelengths. The reason why they are observed in different spatial 
directions is explained by the phenomena analysed in Sections 5 and 6 above -refraction 
and reflection by spherical rain-drops. The angles at which the green and red bows are 
observed are controlled by the vàlues of the scattering amplitude A sca t and the wavelength 
A° in (5.9) causing different wavelength photons to follow diffent paths according to Snell's 
law of refraction, (6.11). 

Now I take a compact disc (CD) and hold it in the direct sunlight, turning it in my 
hand. At certain angles of the face of the CD relative to the direction of the Sun I 
observe vivid flashes of green or red light. As for the rainbow, the photons striking the 
cone receptors in my retina, when I perceive the green or red colour, occur in narrow 
bands of wavelength. However, the physical source of this sensation of colour is yet again 
different. The face of the CD constitutes a diffraction grating. The bright flashes of colour 
correspond to constructive interference of the path amplitudes of photons for particular 
intervals of wavelength. The quantum mechanical description is essentially the same as for 
the photon YDSE described in Section 9, except that, not just two, but many, probability 
amplitudes interfere constructively [74] to give large vàlues of the combined probability 
amplitude Afi in (2.1) and therefore a large signal in the output neurones of the cone 
celi receptors. Where is now the physical origin of the visual perception of the different 
colours? It is not simply the grooves on the surface of the CD. These serve simply to 
define the lengths of the photon paths from their begining, at an atom in the Sun, to their 
end, at a cone celi in my retina. Constructive interference corresponds to a constant phase 
difference between the amplitudes of photon paths reflected between successive grooves 
on the CD. In Feynman's interpretation of quantum mechanics, these phase diffences are 
those between alternative paths of the same photon reflected from different places on the 
surface of the CD in its alternative histories. Constructive interference occurs when, for 
each adjacent pair A,B of paths, 4>b~<Pa = 27rn where n is an integer. The actual physical 
source of this phase difference is, uniquely, the temporal propagator (2. 12), (3. 20) of the 
source atom. Thus, as I move the CD to observe either green or red interference fringes I 
select photons from different sets of source àtoms in the Sun, and am, in fact, observing 
the temporal phase advances of the wavefunctions of these different àtoms (Feynman's 
'stopwatch hand' [1]) as probed by alternative paths followed by each single photon that 
is absorbed in the cone celis of my retina. These phase advances actually occur some 8.3 
minutes before I perceive the fringes! 

The last experiment on the perception of colour to be described is the one featured 
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on the dust-jacket of the first edition of Feynman's book [1], that provided the primary 
motivation for the writing of the present paper. In order to observe this effect -thin film 
interference, under the best conditions, an overcast, rainy, day is desirable. This is for two 
reasons, firstly because wet tarmac road surfaces favour the production of thin films when 
small oil-drops fall upon them from passing vehicles, and secondly because diffuse sunlight, 
multiply scattered and refracted from the wàter vapour droplets of clouds, ensures that 
photons with paths that satisfy conditions for constructive or destructive interference can 
arrive at the eye from all parts of the thin oil-patch. In this case, the perceived colour 
is affected both by refraction, as in the case of the rainbow, and by path differences as 
in the case of the photon YDSE, or the CD experiment. The angular positions at which 
the differently coloured 'Newton's Rings' appear depend on both the thickness of the oil 
film and, as Newton noted [2], on the refractive index of the oil. However, the bàsic 
physical mechanism, in Feynman's interpretation, is the same as for the CD experiment - 
temporal phase advances of the wavefunctions of àtoms in the Sun, several minutes before 
the photons from the decay of those àtoms strike the receptor celis of the retina of the 
observer. 

Indeed, quantum mechanics is essential not only for 'the interpretation of the experi- 
mental material' but to the very physical mechanism of colour vision itself in the eye and 
the brain. The cone celi receptors are quite analagous to the leaves and flowers of the 
geranium i.e. they function essentially as wavelength filters. They do this because of the 
different structures of atòmic energy levels to be found in the three different types of cone 
receptor celis. The world of our every-day experience not only requires quantum mechan- 
ics for its rational interpretation but also for enregistering the raw data from which that 
experience is derived. 
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Appendix A 



The case of a rectangular boundary, of sides Ly, Lz for the refractive medium in 
Fig.2 or Fig.3 is shown in Fig 12a. The Ox axis, along which are situated the excited 
source atom and the photon detector, intersects the YZ plane at Op and has cartesian 
coordinates (Z,Y) relative to the center, Op, of the rectangular boundary. The axes OY 
and OZ are parallel to the sides of the rectangle B^B^B-^B^ whose plane is perpendicular 
to Ox. The typical boundary point, B, has polar coordinates (Ri, 0i) relative to Op The 
corresponding value of r™ ax in (5.2) (see also Fig.2) is given by: 

rr * = ^ + Rl =Xl + R + 0(4) (Al) 

LX\ x^ 

The angles $i, $2, $3 and $4, of the line segments OpBi, OpB2, OpB% and OpB^, 
joining the còrners of the rectangle to Op, are defined, in an anti-clockwise sense, relative 
to the axis OZ,. 

The function r™ ax (di, 0i) in (5.2) then has the following form, where corrections of 
0(R\/x\) are neglected: 

$4 < 01 < 27T + $1 

rr(ll , iz ,Z, W = * 1 + %^f (.42) 



$1 < 01 < $ 2 

rr*(*i, Ly, Y, 00 =x 1+ {L Y,2 .~ 2 Y J (A3) 

2xi sm 0x 



$2 < 01 < $3 
rto,L,,Z M =* l + ££±g (.44) 



$3 < 01 < $4 

rT ax (x 1 , L Y , y, 00 =x l+ {L Y/2+ 2 Y / (Aò) 

2xi sm 0i 

The case of a circular boundary is shown in Fig 12b. By symmetry, it is sufficient to 
consider the displacement, Y, of of the projection, Op, of the Ox axis from the center 
of the circle, Oc, whose plane is perpendicular to Ox. Choosing the OY axis parallel 
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a) b) 

Figure 12: Dependence of the upper límits of the r\ integral in (5.2) on the azimuthal 
angle <pi for diSerent boundary geometries. Op is the projection of the classical rectilinear 
path joining the source and detector on the YZ (transverse) plane. Various geometrical 
parameters used to specify r™ ax are defined for the case of rectangular, a), and circular, 
b), boundaries. 
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to OcOp, and a boundary point, B, with polar coordinates (Ri, <f>i) relative to Op, the 
function r™ ax {di, <f>i) is found, from the geomery of Fig.lOb, to be: 

(JR 2 B -Y 2 cos 2 0! - Y sin 0O 2 
rr x (x 1 ,R B ,Y, ( j )l )=x 1 + K -^ ^ (A6) 

where Rb is the radius of the circular boundary. 

The following contributions of the upper limit of the r\ integral in (5.2) are then found 

Iu = ~ exp[i/tr™ a:c (^,0i)rf0i 
Jo 



2tt 



2 7r r™ ax (d i ,(f) 1 ) . . 2nrV iax (d i , 00 
cos + % sm 



#! (A7) 



where 

— Xi, Ly, Lz,Y, Z rectangular boundary 

di = Xi,Rb,Y circular boundary 

In the typical case that r™ ax ^> À°, a very ràpid phase variation of the integrand in (A7) 
occurs as a function of <f>i, so that the integral Iu over this variable vanishes. 

An possible exception to this is the case Y = in (A6), in which case: 



r? ax (x 1 ,R B A<l>i) = Xi + ïr- (A8) 

lX\ 

which is independent of 0i, so that no cancellations due to ràpid phase variation as a 
function of this variable occur. However, even in this case, both the position of Op must 
be specified, and the boundary radius Rb must be constant, to within a fraction of À° in 
order to avoid the cancellation of Iu due to ràpid phase variation of the sine and cosine 
functions in (A7). 
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Appendix B 

The nested n fold integral T n (L) is defined as: 

l n (L) — / dxi dx 2 ... / dx n (Bl) 

Jxi J X?. Jx—kr 



Consider the incremental change: L — > L + <5L. Substitution in (Bl) gives: 



/ / >:r+ 2" f XA — \ l f X+ 2 f 

I n (L + 5L) = i dx x + / dx 1 / cfe 2 + / 



dx? 



^2 

1X2 Jx +^ / \ Jx 3 Jx + J ï / 

J n (L) + i x + i 2 + ... + ij + ... + i n + i_ + 0(5L 2 ) (52) 



where 



<5L 

ii = — / cfe 2 / tte 3 ... / cfe n = — J n _i(L) (53) 



L+SL - , i 



/ c&ri / dx2-.. / dx,-... / c&r n (54) 

A 2 Jx 3 J x +jr ^-f 

/■s+f r^+í / >í +¥ ^"í 

/ tZxi / dx 2 ... / cb n _i / dx n (55) 

^ £2 X3 J X n X ~ 



2 

Because of the linear dependence of each integral on its upper and lower limits, the 
following replacement may be made in (B4): 



where 



similarly 



where 



f X +2 , f X +2~ , f X +2 ÒL 

.. / dxj-! / dxj... -> ... / — . 

JXj Jx+^ JXj Z 

L 5L 



% = x + + -j 



dx,-_ 2 / — ... -> ... / — — . 

íe,_i Jij Z Jxj — i ^ Z 



= x + — + 



2 

It follows that i 2 , is, ■■■i n are at least of order SL 2 . Similar considerations show that: 

r x +i r x +^ r x +i 5L 

i- — dx\ I dx 2 ... j dx n -i— (56) 

JX2 Jx3 JX n 2 



where 



or 



L 5L 

X n ~ X - - - — 



rx+f rx+f I rx+f rx-f \ SL 

L dXi L dX2 - [L· dXn ~ i + Lk-^ dXn ~ i t 



IX 2 JXz 

Y^n-i(L) + 0(SL 2 ) (B7) 
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2 4 



Since %2i Í3,---i n give no 0(SL) contribution, combining (B2),(B3) and (B7) then gives: 

l n {L + SL) = l n (L) + l n ^ÓL + 0{5L 2 ) (B8) 
or, from the Taylor expansion of l n (L + SL), 

dZn(L) 



hence 



Since 



then 



dL 



= 2 n _ 1 (L) (B9) 



l n (L) = [ L l n ^(L)dL (510) 
Jo 

T\{L) = [ X+2 dx! = L (fill) 

Jx-k 



ML) = 


[ L LdL = 
Jo 


L 2 
2 


UL) = 


/ —dL 
Jo 2 


L 3 
~ 3! 


ML) = 




L 4 
~ 4! 



which is (5.16) of the text. 



UL) = / t — — g?L = - (512) 

jo (n — 1)1 n\ 



79 



Appendix C 

The nested ri,r 2 ,r 3 integral in (4.30) is: 

/•As — T2— r- L í+xi rAs — rz+X2 rAs+X3 

J 3 (r) = / e lKri d ri / e tKr2 dr 2 / e lKr3 dr 3 

Jx\—X2 JX2—X3 Jx3 



Performing the T\ integral: 



IKX1 



—e lKri dn 



L 



As—rz+X2 



gi/í(As— rs) ^k.{t2-X2) 



/■As+X3 . 

dr 2 / e tKrs dr 3 

Jx 3 



Performing the r 2 integral: 

rAs+x 3 



X3 



(Ag-r 3 + x 3 ) cmAs _ (e^ As - e ÍK ^- x ^) 



IK 



(ík)'" 



dr 3 



Finally, performing the r 3 integral: 

3 r / (ik/as) 2 



2o( r = e lKX1 



+ íkAs — 1 e 



íkAs 



(Cl) 



The factor exp[ÍKXi] is cancelled by the factor exp[— íkxi] in (5.30) so that the integrand 
of the Xi, X2 and x 3 integrals in this equation is constant. The integral is then given by 
(5.16) of the text as the r and x integrals factorise. The replacement kAs = A<£> in (Cl) 
yields (5.31) of the text. 

The general expression for the nested n-fold r and x integrals obtained in a similar 
manner to the three-fold case just discussed is, in the notation of (5.33): 



X„(x,r) = 



{0L) n 



nl 



(iA<i>y 



n- 11! 



(n-2)! 



+ iA$ - 1 I e íA<ï> + 1 
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Appendix D 



Multiplying exp[ïA$] by the polynomials in iA<£> at each order in j3L, the expresion 
in large curly brackets of (5.33) may be written as: 



Ve/ 



+ 
+ 
+ 
+ 



1-iPL 

2! 
3! 

r 

4! 

(i(3L) n 



' . (iA$) 2 (iA$) 3 
ïA$ + v - 7 + v - 1 + ... 



2! 
1 1 
ÏÏ ~ 2! 
1 I 



3! 

(iA$) 3 
1 



1 1 
2! ~ 3! 



+ ... 



* A$) ( ~2! + 2! " 3! ] + (ïA$) 



^-^7 + 



(iA$) 4 
(iA$)" 



1 



1 1 1 

+ 777 - 



3! 2!2! 3! 4! 
1 | 



+ (iA$) 
1 



1 1 


1 


2!2! + 3! 


~ 4! 


5 Í 1 - 


1 


V3!2! 


2!3! 


I) 




n! / 





1 

+ TT - 



+ 



+ (iA$) 



+ 



n+l 



1 



(n-l)!2! n! (n + l)! / 



+ 



(Dl) 



+ ... 



Separating out the real and imaginary parts at each order in f3L the gives the expres- 
sion: 



ref 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



l + (3L 
2! 



A$ 



(A$£ . / (A$) 2 

3! + - + ^ 2! 



(A$) 4 



4! 
(A$) 4 



+ ... 



(A$) 3 (A$) f 



3! 



5! 



3! 



+ ... + A$( 



2! 
(A$) 2 
2! 



4! 
(A$) 4 
4! 



(A$) ; 



3! 



(A$) 5 , ^ , A$( (A$) 2 | (AD) 4 



5! 



2! 



4! 



...) + 



(A$) ; 
2! 



(A$ 



(A$) 6 , ^ , A$( (A^) 3 , (A$) 5 



4! 
(A$) 2 



6! 3! 5! 

(A$) 4 + (A$f_ ^ + A$( (A$) 3 



...) + 



(A$) 2 ,(A$) ; 



(A$) ; 
" 3! 
(A$) 4 



+ 



2! 



2! 



4! 



+ .. 



4! 
(A$) 2 



2! 



2! 



+ 



6! 
(A$) 4 



3! 



(Aff 
5! 



+ ...) 



4! 



_...) + í^(A*-í^ + 



+ | j f (A$)5 | (A$)7 



5! 



7! 



... + A$( 



3! 
A<£>) 4 

4! 



3! 



(A$) e 
6! 



(A$) s 
2! 



^ (A$) 3 + (A<£>) 5 



3! 



5! 



^ | (A$) 3 ( (A$) 2 (A$) 4 | j 



3! 



2! 



4! 



+ 



(D2) 



Noting that all the infinite series in A<3> in (D2) are sine or cosine series, possibly with 
sequences of missing low order terms, enables the right side to be written as in the large 
curly brackets of (5.35). 
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Factoring out the terms proportional to S = sin A<£> and C = cos A$ at each order in 
(3L, D2 may then be re-written as: 



ref 



+ 
+ 

+ 

+ 

+ 

+ 
+ 



l + PL[S + i(l-C)] 

m 2 



2! 
3! 

4! 



[C- 1 + ,SA$ + i(5-CA$)] 

-5(1 - + CA$ + i ^C(l - + SA$ - 1 



2! 



2! 



3! 



(A$) ; 



5(1 



(A$) 5 



5! 



_C(A*-í^)!) + S(l-(^ + (AÏ ' )4 ' 



1 - ,S(A<£> - 



3! 

(A$) 3 
3! 
(£>3) 



)-C(l 



2! 
(A$) 2 



+ 



4! 

(A$) 4 
4! j 



The polynomials in A<£> that multiply 5 and C at each order in f3L are the truncated 
sine and cosine series Sj and Cj respectively, defined in (5.36) and (5.37). Substitution of 
these series in (D3) yields then (5.38) and (5.39) of the text. 
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Appendix E 

It is assumed in the discussion of the Michelson interferometer in Section 8 that the 
source atom is at rest. In the case that the source is in motion in the laboratory frame 
the time factors in the complex exponentials of (8.1) and (8.2) must be scaled by the 
factor I/7 where 7 = — (v/c) 2 and v is the source atom velocity, to take into 

account relativistic time dilatation. This is because the correct argument of the space- 
time propagator (see (3.11)) is the proper time r, and At = Aí/7 where t is the laboratory 
time. 



Assuming the Maxwellian distribution: 

dN _ 2 -i 
—— = Cp e p 
ap 



(El) 



for the momentum, p, of the source atom, of mass, M, at absolute temperature, T, where 
p 2 = 2MkT, 
the relation: 



p 2 = 2MkT, enables the distribution of the relativistic parameter 7 to be calculated via 



7 



E 
M 



t/NPTp 2 
M 



, 1 P 2 
l + 2W 2 + 



(E2) 



Modifying the argument of the cosine in (8.7) to take into account the motion of the 
source atom gives 



/n A ,2nd. /n , Kp 2 d. 
cos(2/trf) — > cos( ) ~ cos(2/trf ^2~) 

Performing the average over p using the distribution (El): 



(cos(— )) = Re e 2 * Kd 
7 



v 2 ■ 2 

j e p e m 1 



dp 



p 2 
^2". 



J °°e p dp 

Performing the integrals (see [75]) gives the result: 



(ES) 



(EA) 



<cos(^)> 



Re 



k 2è 



2ind 



7 



p3 



{Eh) 



where 



k = + ík 



_p_y 

M) 



d] 



Since k (j^} d 1, it follows that: 



k 2 ~ p 3 [l 



-ik 



p_ 

M 



d] 



~ p 3 exp 



-IK 



_p_y 

M/ 



(#6) 



Combining (E5) and (E6) gives finally: 



(cos 



7 



-)) 



cos 



P 



(El) 



Thus the only effect of random Maxwellian motion of the source atom is the small relative 
change — (3/4)(p/M) 2 of the phase of the interference term. There is no damping of the 
fringe visibility due to Doppler smearing of the wavelength of light as predicted by the 
classical wave theory. The numerical size of the correction term 3/4(p/M) 2 for a sodium 
atom at NTP is 1.6 x 1(T 12 . 
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Appendix F 



For convenience, units with h — c — 1 are employed in this Appendix. The following 
variables are used to characterise the velocities va and vb of the electron in path A or 
path B: 

' ' /3 = t A -t B (FI) 



ot = 



The velocities are given, in terms of these variables and the path lengths r' + r A and 
r' + r B by the relations: 



[l-a)(r' + r A ) 

(3 



v B = 



(1 -ct)(r' + r B ) 

a/3 



(F2) 



The physical ranges of a and /3 are: 



i <. u; <. . . — ot max 

(r' + r A ) 

< < oo (F4) 



(F3) 



The limit /3 — > oo corresponds to vanishingly small velocity. The vàlues a — 1, (3 — (in 
which case the ratio (1 — a)//3 remains finite) correspond to equal production times for 
any allowed velocity, whereas ot = a max , for any (3, corresponds to equal velocities. The 
equations (F2) give for the velocity difference: 



r'{l-a) 



1 i\ , 1 í r ' B \ 
--!) + -( r A) 

ot r a 



and for the average velocity: 



_ v B + v A r'(l - a) 
v = = 



2(3 



ot 



r ot 



(F5) 



(F6) 



Combining (F5) and (F6): 



vb - v A 



2v 



G + V + HT + rA) 
{r' + r A )(^^- 



(r 1 + r) 
25 v + 0[(Ar) 2 } 



+ 0[(Ar) 2 ] 

(F7) 



To first order in Ar 



= V 

From the relativistic relation: 

it follows that S v = (m 2 /E 3 )5, 

Pa = P 



■S v , v B = v + S v (F8) 
v= P_ = P 

E yj(m 2 +p 2 ) 

so that 

■à P , Pb =P + à p (F9) 
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where 



p(l + (|) 2 )(r- + r A )(^-l) 
°p m 2°v 2(r' + r) 



(FIO) 



Substituting now (F9) into the formula (9.15) for the phase difference, retaining only first 
order terms in 5 p , and using (FIO), gives: 



+ 2 r' + rH 



P 



Ar / p l\ . a r 



-1^ 



P 



(Fil) 



Substituting a = a max in (Fil), corresponding to va — vb, gives 

m 2 Ar 



<t>B ~ <Pj 



P 



(equal velocities) 



(F12) 



while setting a = 1 or ía = íb gives: 

0b — 0A = pAr (equal production times) (-^13) 

For non-relativistic partides for which p <^ m the equal velocity relation (F12) gives an 
effective de Broglie wavelength much shorter than (9.19) that is applicable in the equal 
production time case, (F13), and which is in good agreement with experiment 
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